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Preface 


This  report  is  the  first  of  three  volumes  which  are  the  product  of  an  analytical 
and  experimental  investigation  into  the  effect  of  a  hydraulically  supported, 
tilting-pad,  Journal  bearing  on  the  attenuation  of  noise  originating  from  rotor 
unbalance. 

The  study  is  a  part  of  a  larger  program  to  reduce  structure -borne  noise 
originating  in  high-speed  rotating  equipment  onboard  Navy  vessels.  The 
introduction  of  additional  flexibility  at  the  Journal  bearings  was  proposed  to 
reduce  noise  transmission,  and  in  the  case  under  study  the  flexibility  is 
provided  by  hydraulic  pistons  and  accumulators.  In  order  to  evaluate  and 
optimize  the  additional  flexibility,  this  is  a  comprehensive  investigation  of 
the  dynamics  of  a  rotor-bearing  system  and  how  the  system  behavior  is  affected 
by  the  dynamic  properties  of  both  the  hydrodynamic  oil  film  and  the  flexible 
supports. 

The  second  volume  of  this  study  will  be  concerned  with  the  force  transmission  of 
(l)  two-mass  rotor  bearing  system  and  (2)  uniform  rotor  bearing  system.  The 
third  volume  will  include  (l)  general  rotor  analysis  and  (2)  experimental 
correlation. 

Mechanical  Technology  Incorporated  was  primarily  responsible  for  the  analytical 
portion  of  t  he  investigation  while  Westinghouse  Electric  Corp.  designed  and 
conducted  the  experimental  test. 
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ABSTRACT 

A  method  for  calculating  spring  and  damping  coefficients  for  the 
tilting  pad  journal  bearing  is  presented.  The  analysis  includes  the 
effect  of  pad  Inertia.  Numerical  results  are  given  in  form  of  de¬ 
sign  curves  for  the  centrally  pivoted  4  shoe,  5  shoe,  6  shoe  and  12  shoe 
journal  bearing.  A  comparison  with  test  results  is  included. 


INTRODUCTION 

This  report  is  part  of  a  study  of  the  force  transmitted  from  a  rotating 

* 

shaft  to  its  bearing  pedestals  where  the  force  is  generated  by  the  presence 
of  an  unbalance  in  the  rotor.  The  shaft  is  supported  in  fluid  film  bear¬ 
ings  and  since  the  fluid  film  possesses  both  flexibility  and  damping  the  re¬ 
sulting  rotor  motion,  and  therefore  the  transmitted  force,  are  directly  in¬ 
fluenced  by  the  properties  of  the  bearing  film.  Thus  it  is  necessary  to 
determine  the  bearing  stiffness  and  damping,  expressed  in  terms  of  spring 
and  damping  coefficients,  in  order  to  calculate  the  transmitted  force.  The 
bearing  spring  and  damping  coefficients  depend  on  the  operating  conditions 
and  the  bearing  configuration  and  have  been  determined  for  several  bearing 
types  in  a  preceding  investigation,  see  Ref.  2.  In  connection  with  the 
present  conttiinuation  of  the  investigation  it  has  become  necessary  to  examine 
yet  another  bearing  type,  namely  the  tilting  pad  bearing. 

References  1  and  2  show  that  although  proper  bearing  design  procedures  can 
effect  a  reduction  in  the  transmitted  force,  bearing  instability  imposes 
an  upper  limit  on  the  attainable  force  attenuation.  Hence,  to  further  im¬ 
prove  the  attenuation  other  methods  must  be  employed.  Two  possibilities  seem 
open:  a)  the  use  of  an  inherently  stable  bearing  type  such  ae  the  tilting 
pad  bearing,  or  b)  by  vibration  isolation  of  the  bearing  housing.  In  order  to 
achieve  a  substantial  force  attenuation  the  latter  possibility  must  be  chosen. 
However,  this  suggests  a  very  soft  bearing  support  which  in  turn  implies  a 
potential  risk  for  bearing  instability.  Hence,  from  both  points  of  view  the 
tilting  pad  journal  bearing  offers  advantages.  For  these  reasons  »n  analysis 
of  this  bearing  has  been  undertaken  and  it  is  the  purpose  of  the  present  report 
to  give  the  results  of  the  investigation  in  form  of  design  curves  for  the 
spring  and  damping  coefficients.  A  second  report  will  deal  with  the  force 
attenuation  due  to  the  flexible  bearing  support. 


RESULTS 


Numerical  calculations  have  been  performed  for  four  bearing  geometries: 
the  4-shoe,  the  5-shoe,  the  6 -shoe,  and  the  12-shoe  tilting  pad  journal 
bearing.  All  four  configurations  have  centrally  pivoted  pads.  The  re¬ 
sults  are  given  in  form  of  graphs,  Figs.  1  to  19,  where  the  spring  and 
damping  coefficients  are  given  as  functions  of  the  bearing  Sommer f eld 
number.  The  following  factors  are  investigated: 

a)  bearing  length-to-diameter  ratio 

b)  the  direction  of  the  static  load 

c)  the  inertia  of  the  shoes 

d)  the  preload  of  the  bearing  (for  vertical  rotors) 

The  results  are  given  in  dimensionless  form: 

CK  CK 

dimensionless  spring  coefficient:  xx  ,  yy 

W  W 

CcdC  QdC 

dimensionless  damping  coefficient:  xx  ,  vv 

W  W 

Sommerfeld  Number:  S  =  (77)  ^ 

W  C 


where: 

C 

Pad  radius  of  curvature  minus  journal  radius,  inch 

W 

Bearing  reaction,  lbs. 

D 

-  Journal  diameter,  inch 

R 

-  Journal  radius,  inch 

L 

Bearing  length,  inch 

2 

-  Lubricant  viscosity,  lbs. sec/in. 

N 

-  Rotational  speed,  RPS 

to 

2jtN ,  angular  speed,  rad/sec. 

K  ,K 
xx  yy 

-  Spring  coefficients,  lbs/in. 

C  ,C 
xx  yy 

Damping  coefficients,  lbs. sec/in. 

Figures  1  to  12  are  intended  as  design  curves  for  horizontal  rotors. 

They  cover  the  three  most  commonly  used  bearing  geometries  for  three 
values  of  the  L/D-ratio.  The  first  9  graphs  have  the  static  load  direc¬ 
tion  between  the  bottom  pads  whereas  for  the  next  3  graphs  the  load  vector 
passes  through  the  pivot  of  the  bottom  pad.  Since  the  bearing  geometries 
are  symmetric  about  the  vertical  load  line  and  shoe  inertia  is  neglected 
the  cross -coup ling  spring  and  damp  coefficients  vanish.  Furthermore,  for 
the  4-shoe  bearing  there  is  also  symmetry  around  a  horizontal  line  such  that 
coefficients  in  the  horizontal  and  the  vertical  direction  are  identical,  i.e. 

K  =  K  and  C  «■  C  This  does  not  hold  for  the  5-shoe  and  the  6-shoe 

xx  yy  xx  yy 

bearing. 


On  each  graph  is  shown  a  curve  labeled  "Critical  Hass".  It  gives  the  value 
of  the  pad  inertia  necessary  to  incur  resonance  of  the  pad  motion,  resonance 
being  defined  by  a  90  degree  phase  angle  between  journal  and  pad  motion  (or 
in  other  words,  at  resonance  the  phase  angle  between  a  zero-inertia  pad  and 
a  pad  with  a  "critical  mass"  is  90  degrees).  Hence,  for  a  particular  bearing 
design  the  value  of  the  dimensionless  "pad  mass"  can  be  calculated: 


where: 


dimensionless  mass 


CWM 

(*Mt>2)  2 


X  lbs • sec 

M  =  o2  equivalent  pad  mass,  — r1 - 

Rp  in.  2 

I  =  Transverse  mass  moment  of  inertia  of  pad,  lbs. in. sec 
Rp  =  Pad  radius,  inch. 

W  =  Bearing  reaction,  lbs. 

D  =  Journal  diameter,  inch 

R  =  Journal  radius,  inch 

L  =  Bearing  length,  inch 

C  =  Pad  clearance,  inch 

2 

ia  =  Lubricant  viscosity,  lbs.sec/in. 


Entering  the  appropriate  graph  with  this  value  the  corresponding  "critical" 
Sommerfeld  number,  and,  therefore,  the  resonant  speed,  can  be  determined. 


Figures  13  to  16  illustrate  the  effect  of  pad  inertia  on  the  spring  and 
damping  coefficients  for  a  particular  bearing  geometry,  namely  the  4-shoe 
bearing,  L/D  ■  .75,  with  the  load  passing  between  the  two  bottom  pads. 

The  calculations  are  performed  for  six\alueB  of  the  dimensionless  pad  mass. 
Since  the  pads  are  given  inertia  the  cross -coupling  terms  no  longer  vanish 
as  for  Figs.  1  to  12. 


Figures  17  to  19  apply  to  a  vertical  rotor  with  zero  bearing  eccentricity 
ratio.  In  order  to  have  stiffness  the  bearing  must  be  preloaded.  The  pre¬ 
load  is  defined  by: 

Preload  =  (1  -  C^) 

where: 

C  »  difference  between  radius  of  curvature  of 
pad  and  journal  radius,  inch. 

C1  ■  difference  between  radius  of  pivot  point 
circle  and  journal  radius,  inch. 

Hence,  C ^  =  1  corresponds  to  no  preload  whereas  »  0  means  that  the 
journal  touches  the  pad. 


Since  the  bearing  reaction  W  is  zero  for  a  vertical  rotor  the  coefficients 
must  be  made  dimensionless  in  a  form  different  from  the  horizontal  rotor 
results.  The  following  form  is  chosen: 


vertical 

rotor 


dimensionless  spring  coefficient 


dimensionless  damping  coefficient 


dimensionless  critical  mass 


HNL(*)  3 
coC 

XX 

hnl(*)3 


M 


j  .  R\  3 


where  the  symbols  are  defined  above. 
K  »  K  and  C  =  C  . 


xx  yy  xx  yy 


Note,  that  for  a  vertical  rotor 


All  the  given  results  are  based  on  partial  arc  fluid  film  force  deri¬ 
vative  computed  by  finite  difference  calculations  on  a  computer  using 
program  PN0091,  Appendix  C.  Fluid  film  rupture  is  included.  Each 
partial  arc  (in  total  10)  is  calculated  for  8  eccentricity  ratios:  e  **.01, 
. 1,  .2,  .35,  .5,  .65,  .8  and  .95.  A  summary  of  the  derivatives  is  given 
in  Table  1. 


DISCUSSION  OF  RESULTS 


A  comparison  between  calculated  and  measured  results  is  shown  in  Figs. 

20  and  21.  The  test  results  are  taken  from  Ref.  4.  Even  if  Ref.  4 
gives  experimental  data  for  values  of  the  Sommerfeld  number  ranging  from 
0  to  16,  Ref.  5  indicates  that  the  actual  measurements  are  limited  to  a 
range  from  .15  to  1.2  and  that  outside  this  range  extrapolation  has  been 
used.  Based  on  these  considerations  the  agreement  seems  reasonably  good 
within  the  normal  range  of  bearing  operation.  It  should  also  be  noted  that 
the  test  measurements  are  obtained  with  vibration  amplitudes  of  up  to  50 
percent  of  the  minimum  film  thickness  whereas  the  theoretical  calculations 
assume  very  small  amplitudes. 

The  most  unusual  aspect  of  the  theoretical  results  is  the  sudden  loss  of 
bearing  stiffness  when  approaching  pad  resonance.  However,  the  effect 
may  be  less  drastic  than  shown.  Since  at  resonance  the  "critical  mass" 
is  zero  the  effect  of  pad  inertia  cannot  be  neglected.  Then  Figs.  13  to  16 
show  that  the  cross-coupling  terms  become  important.  Hence,  the  overall  co¬ 
efficient  may  not  be  zero  (see  Eq.  (10),  Ref.  1).  Pad  resonance  occurs  in 
general  at  high  Sommerfeld  numbers  where  the  eccentricity  ratio  is  small  and 
thus,  the  out-of-phase  between  pad  and  journal  motion  may  not  be  too  serious. 
In  addition  the  fluid  film  damping  is  large. 

Finally,  it  should  be  noted  that  the  fluid  film  force  calculations  include 
film  rupture.  When  the  bearing  is  flooded  with  oil  the  film  can  sustain 
pressures  below  the  ambient  and  rupture  does  not  take  place.  Instead  the 
film  may  cavitate.  Therefore,  if  the  present  results  are  applied  to  the 
calculation  of  a  flooded  bearing  there  may  be  minor  errors  at  low  Sommer¬ 
feld  numbers. 


NUMERICAL  EXAMPLE 

Let  a  bearing  have  a  diameter  D  =  6  inch  and  length  L  *  3  inch,  i.  e., 

L/D  ■  .5.  The  radial  pad  clearance  is  C  ■  .0045  inch  and  the  pad  center 

coincides  with  the  bearing  center  so  that  there  is  no  preload.  The 

speed  is  3600  RPM,  i.e.,  N  =  60  RPS,  and  the  corresponding  oil  viscosity 

*6  2 

is  14  centipoise,  i.e.,  p  =  2.0  •  10  lbs. sec/in  .  With  a  bearing  re¬ 
action  W  =  4,000  lbs.  the  Sommerfeld  number  becomes  S  ■  .244.  When  the 
load  is  between  the  pads  and  pad  inertia  is  neglected  the  dynamic  co¬ 
efficients  for  the  4-shoe  bearing  are  obtained  from  Fig.  7. 

1.15  K  -  1.02  •  106  lbs/inch 
xx 

5.4  C  =  12,700  lbs.sec/inch 

xx 


CK 

xx 

W 

Co£ 

xx 


W 


-3 

The  dimensionless  critical  mass  is  5.25  •  10  which  gives  a  critical 

2 

pad  mass  moment  of  inertia  I  ■  .69  lbs.  in. sec  .  This  corresponds  approxi¬ 
mately  to  a  shoe  weight  of  180  lbs.  which  is  much  larger  than  the  actual 
weight.  Hence,  there  is  no  danger  of  resonance.  For  a  6-shoe  bearing  the 
results  are  obtained  from  Fig.  2  as: 


K  =  4.3  •  10  lbs /inch 
xx 


K  *  1.4  •  10  lbs/inch 

yy 


C  =  10,400  lbs.sec/inch  C  «*  3,400  lbs.sec/inch 
xx  yy 


critical  pad  moment  of  inertia  =  18.5  lbs. in. sec 
corresponding  to  a  shoe  weight  of  12,500  lbs. 


DISCUSSION  OF  ANALYSIS 


The  theoretical  equation  governing  the  fluid  film  behavior  is  the  well- 
known  Reynolds  equation.  Due  to  its  complexity  it  is  normally  solved 
numerically  on  a  computer  and  the  resulting  force  is  a  non-linear  function 
of  the  eccentricity  ratio  and  the  journal  center  velocity.  For  an  exact 
solution  Reynolds  equation  should  be  solved  simultaneously  with  the  equa¬ 
tion  of  motion  for  the  rotor,  leading  to  quite  involved  calculations. 

For  practical  purposes  it  is  convenient  instead  to  assume  that  the  rotor 
amplitude  is  sufficiently  small  to  allow  replacing  the  fluid  film  forces 
by  their  gradients  around  the  steady  state  operating  eccentricity.  Thus, 
the  forces  become  proportional  to  the  vibratory  amplitude  and  velocity, 
the  coefficients  for  proportionality  being  denoted  spring  and  damping 
coefficients.  Introducing  an  x-y-coordinate  system  the  fluid  film  forces 
can  be  written  as: 


F  =  -K  x  -  C  x 

-  K 

y  -  C  ' 

X  XX  .  XX 

xy 

xy 

F  =  -K  x  -  C  x 

-  K 

y  -  C 

y  yx  yx 

yy 

}  yy 

such  that  there  are  A  spring  coefficients  and  4  damping  coefficients. 

For  a  fixed  pad  the  spring  and  damping  coefficients  can  be  computed  from 
the  gradients  of  the  fluid  film  force  as  shown  in  the  Analysis  Section  or 
in  Refs.  1,2,3  and  6.  The  gradients  are  obtained  by  numerical  differen¬ 
tiation  of  computer  results. 

The  analysis  of  the  tilting  pad  bearing  assumes  the  spring  and  damping 
coefficients  for  the  fixed  pad  to  be  known.  Thus,  for  an  arbitrary  rotor 
motion  it  becomes  possible  to  establish  the  equation  of  motion  for  the  pad 
itself,  including  the  inertia  of  the  pad.  Coupling  this  equation  with  the 
above  equations  for  the  fluid  film  force  yields  the  spring  and  damping  co¬ 
efficients  for  the  tilting  pad.  A  summation  over  all  pads  results  in  the 
combined  spring  and  damping  coefficients  for  the  complete  tilting  pad  journal 
bearing. 
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ANALYSIS 

Steady  State  Egullbrlum 

Referring  to  Fig.  22,  let  the  center  of  the  tilting  pad  bearing  be  Og. 
The  bearing  is  made  up  of  a  number  of  tilting  pads,  the  arbitrary  pad 
having  the  pivot  point  P  located  in  an  angle  from  the  vertical  load 
line.  The  pad  is  free  to  tilt  around  the  pivot  point  which  for  conven¬ 
ience  is  assumed  to  be  located  on  the  surface  of  the  pad.  The  steady 
state  position  of  the  journal  center  is  0^  which  is  the  origin  of  two 
fixed  coordinate  systems:  the  x-y-system  (x-axis  vertical  downward, 
y-axis  horizontal)  and  the  ^-o^-system  (the  £-axis  passing  through  the 
pivot  point).  The  location  of  03  with  respect  to  the  bearing  center 
is  given  by  the  eccentricity  Og03  =60»Cce  and  the  attitude  angle  (f)t . 
The  steady  state  position  of  the  pad  center  is  designated  0h  such  that 
the  journal  center  eccentricity  with  respect  to  the  pad  is  determined 
by  0h  Oj  a  e  »  C  €  .  The  corresponding  attitude  angle  is  (f>  .  The 
journal  radius  is  R,  the  radius  of  the  pad  is  0hP  =  R  +  C  and  the 
radius  of  the  circle  passing  through  all  pivot  points  with  center  in  Og 
is  0BP  *  R  +  C'. 


The  point  0no  is  the  pad  center  with  no  tilting.  Hence,  0h#  0„  is 
circular  arc,  or  for  small  motions,  a  line  perpendicular  to  On«P. 
Projecting  0-j  on  0h,P  yields: 


ecosf-  I-  c  ~ iecosiy-f0) 


This  equation  contains  three  unknowns , namely  6,  and  (po  since  is 
the  independent  variable.  The  second  equation  derives  from  the  requirement 
that  the  force  on  the  pad  passes  through  the  pivot  point  which  establishes 
a  relationship  between  €  and  (j?  (i.e.  the  journal  center  locus  with  re¬ 

spect  to  the  pad) .  The  third  relationship  is  the  requirement  that  the  total 
horizontal  force  component  (in  the  y-direction) ,  summed  over  all  pads,  is 


Z  Fs.'hTJ/  =0 

aU  1 


(2) 
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Eq.  (2)  is  usually  used  to  determine  (y0  by  trial-and-error  as 

ft  ' 

follows:  for  a  particular  case  €0/  £  and  Tp  are  known  in  Eq.(l). 

For  several  assumed  values  of  (j)t  calculate  C(o$^P  for  each  pad, 
determine  the  pad  forces  F  from  available  pad  data  and  plot  ZFsinf 
as  a  function  of  (j>t  .  The  zero-point  determines  the  desired  value 
of  (f6 . 


This  procedure  is  at  best  tedious.  It  becomes  very  complicated  if 
the  pivot  point  is  not  in  the  center  of  the  pad  in  which  case  the 
pad  force  F  can  be  a  multivalued  function  of  ECo5(f).  However,  when  the 
pivot  points  are  located  symmetrically  with  respect  to  the  vertical 
load  line  through  the  bearing  center  Og  then  (|P0  -  0  .  A  further 
simplification  arises  when  the  pivot  point  is  in  the  center  of  the 
pad  since  the  pad  force  F  then  is  uniquely  determined  by  £(oS^. 

These  two  conditions  apply  to  all  the  numerical  calculations  in  the 
present  report  but  is  not  a  necessary  assumption  for  the  analysis  to  be 
valid.  The  analysis  requires  only  that  £  and  (p  are  known  for  each  pad. 
How  these  values  have  been  arrived  at  is  immaterial  in  so  far  as  the 
analysis  is  concerned. 


Fixed  Pad  Coefficients 


In  order  to  determine  the  spring  and  damping  coefficients  for  the  complete 
tilting  pad  bearing  it  is  necessary  to  know  the  forces  and  their  deri¬ 
vatives  for  each  pad  as  if  the  pad  was  fixed.  Under  steady  state  condi¬ 
tions  the  journal  center  has  the  eccentricity  ratio  £  and  the  attitude 
angle  (p  with  respect  to  the  pad  center.  The  fluid  film  pad  force  has 
the  components  Fy.  and  F^  and  under  steady  state  conditions  the  re¬ 
sultant  force  passes  through  the  pivot  point,  i.e.  Fy  *-F  and  F^  «  0, 
see  Fig.  22.  Thus  F  denotes  the  load  on  the  pad.  Usually  the  force  is 
resolved  along  the  radial  and  the  tangential  directions  with  the  components 
Ff  and  Ft,  respectively,  see  Fig.  22.  Hence: 


! 


I-F 

0 


Cosf  Sfh<p 
5ih(f>  -CoS<P 


Ft 


(3) 


For  an  infinitessimal  small  motion  around  the  steady  state  position 
the  dynamic  forces  become  from  Eq.  (3): 


(4) 


0)5  f 

S!n(f> 


Sih(jP 

‘  < 

dFf.  +  d<f> 

- cos(j) 

dFt  -  Fr  df 

The  infinitessimal  dynamic  motion  of  the  journal  center  is  described 
by  the  coordinates  (  ): 

f  =  d(ecos(p) 


or 


(5) 


'dt  ] 

: 

COS^P  Sirtf 

edf 

- sih(f  C0S(f 

*1. 

The  velocities  transform  similiarly: 

(6) 


de 

ii 

<os(f 

sih<y> 

>  < 

• 

! 

edf 

-Slhf 

Cos(f>, 

» 

V 

The  dynamic  force  components  dFr  and  dFt  may  be  expressed  in  terms  of 
the  dynamic  amplitudes.  From  Reynolds  equation  it  can  be  shown  that  the 
fluid  film  force  F  can  be  written  (Refs.  1,2): 

(7)  F* 
where: 

(8) 

(9) 

(10) 

Therefore: 

(11)  «IF-k.{(l-2i)[$Je+$£<lf  +a7S?ft-z4)>l<('4l)]'5££‘'^ 


T. 


c  -  yu N/DL 

V  F  VC 


( Sotnncrfclcl  Nu*>ktr) 


Now: 


4  • 

because  at  the  equilibrium  position  £=(jP-0.  Hence,  Eq.(ll) 
reduces  to: 


Eq. (12)  applies  to  both  Ff  and  Ft>  Let  the  corresponding  dimensionless 
forces  be  denoted  f  and  as  defined  by  Eq.  (7).  Thus,  by  substitution 
of  Eq.  (12)  into  Eq.  (4): 


fl|f^  COSf  sih<p 

(13)  J 

JR]  5  inf  -COSf 


d£  /ift  +  it) 

/ift  -it) 
d€  Udf  t) 


Define  the  fixed  pads  spring  and  damping  coefficients  by: 

*  -^fjr f  “  ^rr i  ~  ^r>  "1?  ~  ^ 

(14) 

To  determine  the  8  coefficients,  substitute  Eq.  (5)  and  (6)  into  Eq.(13) 
and  collect  the  terms  in  accordance  with  Eq.(14)  to  get: 

(15)  Kjj-fWf^Coi'ip-  E^p 

(16)  as'f  +  Hi  Mfsinlf  -  ^  S/n<p] 

(17)  «*’<r  $  + 

(18)  [j|j  sin 'If  +  jfe  C„S(}Si»<?  +  ¥  (osf  1 

(19)  ViM-f  «*P  -  tf>  ‘"Wi*  +  «  £  Si"fl 

(20)  coC1s=£A<o  ['Hi  tos’f  +  ^  (os/psinf  -  &  snip] 
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( 21 )  K,,  - i  ><o [- jjf  ta'f  +  $  Sin’f  +  -  li* ) »sf  sinf  ~  £  as f  ] 

(22)  «CV,=£A*»[  sm*<p  -  H c«fs.nf  +  ^  »«f] 

where: 

(23)  (s 

(24)  ^  *  0 

(25)  C^40=  c  PSp 


and  all  forces  and  derivatives  are  calculed  for  the  given  steady  state 
position,  defined  by  e. 

Tilting  Pad  Coefficients 

Referring  to  Fig.  22,  0n  is  the  steady  state  position  of  the  pad  center. 

Under  dynamic  load  the  pad  center  oscillates  around  0^  with  the  amplitude 
<yjf  such  that  ^p/(?p  represents  the  dynamic  tilting  angle  of  the  pad 
(Rp  is  the  distance  from  the  actual  pad  pivot  to  the  pad  center).  The 
moment  on  the  pad  from  the  fluid  film  pressure  is  “■  I?p  dF*f  .  Therefore, 
if  the  mass  moment  of  inertia  of  the  pad  is  I  the  equation  of  motion  becomes: 

ifr-w 

or 

(26)  MTrw_<<S 

where: 

M'i 


(27) 
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Thus ,  under  dynamic  conditions  should  be  replaced  by  (  ^-^p)  In 
Eq.  (14).  In  order  to  eliminate  substitute  the  expression  for  dF^ 

into  Eq.  (26): 

(28)  +Cirf  + 

Let  the  dynamic  motion  of  the  journal  center  around  the  steady  state 
position  be  harmonic: 

(29)  ($,'>/)  e‘wt 


Hence,  solve  Eq.  (28): 


<3°)  r'Vf' 


where: 


(31)  P  ‘ 

_  _ _ 

(32)  ftC^-Mw1)**  (aO*)1 


Thus,  replacing  ^  by  Eq.  (14)  becomes: 

(33) 

^  ~  itoGjyly 

where  Ky^  )  etc,  are  the  spring  and  damping  coefficients  for  the 

tilting  pad  and  given  by: 

(34)  Kfj= 

(35)  ^  Cjj  »  <X>C^-(p kj^+^0)C^) £oC^  +(q  k^-ptuCy*,) 

(36)  Kf1*-M<o'l(pk'n+^C^) 
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(37) 

(38) 

•<V- 

-Mw’fpk^j 

(39) 

(oc;r 

Modify 

-pw(V 

(40) 

k*- 

-  Mcoz(pK^ 

■H<o*0+pMw’) 

(41) 

CO 

Mcofyk^ 

-pcoC^)» 

If  the  pad  has  no  inertia,  i.e.  M=>0  t  only  and  CoCf£.  remain 

as  would  be  expected. 


Bearing  Spring  and  Damping  Coefficients 


Having  determined  the  spring  and  damping  coefficients  for  the  individual 
tilting  pads  it  remains  to  combine  them  into  the  overall  bearing  co¬ 
efficients.  The  coordinate  system  for  the  bearing  is  the  x-y-system, 
see  Fig.  22,  with  the  coordinate  transformation: 


T 

'  COS^ 

smTp 

(42) 

*1 

~ 

-sihf 

cosy 

CoS  Ip 

"Sinlp 

(43)  < 

■-  —  < 

Sihlp 

cosy 

where  dFx  and  dF^  are  the  dynamic  forces  measured  in  the  x-y-system. 
The  bearing  spring  and  damping  coefficients  are  defined  by: 

dfx~~  k’y,,  X  —  C**  X  ~  Kxj  ~  CXI}  tj 

(44) 

dfij—  y  —  Cijx x  *“  tj  —  <j 

Thus,  substituting  Eq.  (42)  and  (43)  into  Eq.  (33)  and  grouping  terms 
in  accordance  with  Eq.  (44)  yields: 

(45)  If,,-  cos'll'  4  K^S.VD'-df^+^lfosy#^ 


(46) 


loC*,  -  UjC‘n  cos'y+eoC^  stff-feoCfr  +  eoC^ )  cos^sinf 

(47 )  Kj  c  -  k^j  SlV  V  +  ( Kft  -  k^  )  (0%  y 

(48)  bC^cofy-coC^  sin'W +(u)Cff-coC'^) cosysihy 

(49)  k,*  -  k^y  coj'^p  -k^  siV^  +  (KfF  - k^)  <<*y  si*y 

(50)  coC^s  cvtijj  Co?\y  -toC^  Sin'f  4  (WCj|  -(oC'^)cosy  sih  y 

(5D  k„*  cosly  sivy-*- (k^+k^) ^y  s^y 

(52)  CoCijy*  CoC^  CoSXy~l-  U)  Sinry  +  (eoCy^  + Codjj)  Cosy  Sihiy 

A  summation  over  all  the  pads  making  up  the  bearing  gives  the  bearing 
spring  and  damping  coefficients.  If  the  pads  have  no  inertia  the 
equations  simplify  to: 

K„  =  %  cos'V  coC%K  -  coC'Jr  co?y 

(54)  kyi^  -  kijjf  -  k yy.  Coi y  Sin'll  4jC)(^  -  toC^  -  CoCyy  Cosy  iihl^ 

kijij  ~  kyy  $ inry  CoCijij  ~  Cyy  Sihxy 

Thus,  for  symmetry  around  the  x-axis  and  no  pad  inertia  the  cross¬ 
coupling  terms  disappear. 


Pad  Motion 


The  pad  motion  is  given  by  Eq.(30)  which  can  also  be  written: 


(55) 


°7r 


_  (ki)j4<U;C^y)f  4(k>«i><|-«-  l  to  )  *1 


M  0ji*"  4  i  CO  C*j‘ 


Let  for  M  =  0,  i.e.  for  no  pad  inertia: 


4  i'to  C 


>=7+S,  f 


(56) 
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Eq.  (58)  gives  the  amplitude  ratio,  i.e.  the  magnification  factor,  and 
Eq.  (59)  gives  the  phase  angle  lag  with  respect  to  a  inertialess  pad. 

Thus  the  two  equations  indicate  how  well  the  pad  follows  the  shaft  motion. 
The  phase  angle  becomes  90°  when: 

(60)  Me,;t  O)2  = 

Hence,  if  the  pad  mass  satisfies  Eq.  (60)  there  will  be  a  resonance  of  the 
pad  motion. 

It  is  convenient  to  use  Eq.(60)  to  establish  a  value  for  M,  designated  the 
critical  mass.  In  dimensionless  form  Eq.  (60)  may  be  written: 

CWM„;(  |  fift*)**  ( 

<61)  [mDlOT  '  W’S  Zga 

where  S  is  the  bearing  Sommerfeld  number. 


CONCLUSIONS 


1.  An  analytical  method  has  been  established  for  calculating  the 
spring  and  damping  coefficients  of  the  tilting  pad  journal  bearing. 
The  coefficients  can  be  used  directly  in  computing  the  critical 
speed,  the  response  and  the  transmitted  force  of  a  rotor. 

2.  Numerical  results  for  several  bearing  configurations  have  been 
obtained  and  are  presented  in  form  of  design  curves.  A  comparison 
with  test  results  shows  fair  agreement. 

3.  The  results  indicate  a  sudden  reduction  in  stiffness  when  approach¬ 
ing  resonance  of  the  pad  motion.  The  implications  of  this  behavior 
have  not  been  assessed. 


RECOMMENDATIONS 


1.  Since  the  force  attentuation  obtainable  with  a  tilting  pad  journal 
bearing  is  not  necessarily  limited  by  oil  whip,  a  study  should  be 
undertaken  on  how  to  optimize  the  attenuation  by  a  proper  choice 
of  bearing  dimensions,  notably  the  clearance. 

2.  The  present  analysis  considers  the  motion  of  the  shoes  around  an  axial 
axis,  i.e.  the  "rolling"  of  the  shoes,  a  more  complete  investigation 
should  incorporate  both  "pitch"  and  "yaw"  of  the  shoe. 

3.  Although  a  tilting  pad  journal  bearing  is  considered  "inherently  stable" 
this  holds  true  only  in  the  idealized  case  where  the  shoes  have  no  Inertia 
and  the  pivots  are  frictionless.  Hence,  the  stability  limit  of  a  tilting 
pad  bearing  should  be  studied,  especially  the  stability  of  the  shoe  motion. 
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APPENDIX  A 


COMPUTER  PROGRAM  PN0078:  SPRING  AND  DAMPING 
COEFFICIENTS  FOR  THE  TILTING  PAD  JOURNAL 
BEARING 


APPENDIX  A 

Computer  Program  PN0078:  Spring  and  Damping  Coefficients  for  the 
Tilting  Pad  Journal  Bearing _ 

The  program  calculates  the  spring  and  damping  coefficients  for  the 
tilting  pad  journal  bearing  based  on  the  analysis  given  in  this  report. 
The  input  to  the  program  consists  of  the  partial  arc  bearing  forces 
and  force  derivatives,  the  position  of  the  pads  and  the  steady  state 
bearing  eccentricity  ratio.  The  forces  and  their  derivatives  can  be 
obtained  from  PN0091,  Appendix  C,  and  d4ta  reduced  by  PN0131,  Appendix  B. 
The  output  from  the  latter  program  serves  directly  as  input  to  PN0078. 

Input  Data 

The  program  is  written  for  the  IBM  1620  computer,  40K  memory  storage, 
with  input  and  output  on  punched  cards. 

Card  1  and  2 

Descriptive  te*t,  Col.  2  to  52 

Card  3 
7  •  (15) 

Word  1  (NP)  gives  the  number  of  pads.  NP  ^  12 

Word  2  (NEB)  gives  the  number  of  bearing  eccentricity  ratios  for  which 
calculations  are  performed.  There  is  no  upper  limit  on  NEB. 

Word  3  (NM)  gives  the  number  of  calculations  performed  for  each  bearing 
eccentricity  ratio  with  different  values  of  the  pad  inertia.  The  maxi¬ 
mum  value  of  NM  is  12.  If  the  pad  inertia  is  neglected,  set  NM  *  1. 

Word  4  (NF)  controls  the  meaning  of  the  input  items  f^  and  ft  in  the  pad 
data  (fr  =  radial  force  component,  f  =  tangential  force  component). 

If  NF  =  0,  then  the  input  value  for  f  is  the  pad  Sommerfeld  number  and 
the  input  value  for  ffc  is  the  horizontal  pad  force  f ^  (usually  f^  *  0). 
If  NF  +  0,  then  <«r)inpJlt«nd  <Vlnput  *  V 
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Word  5  (MC)  controls  the  way  in  which  the  pad  inertia  data  is  given 
as  input.  If  MC=  -1,  each  pad  will  be  given  its  own  inertia.  If  MC  ■  0, 
each  pad  has  zero  inertia.  If  MC  =  +1,  all  the  pads  have  the  same  in¬ 
ertia. 

Word  6  (MS)  controls  the  form  in  which  the  pad  inertia  is  made  dimensionless 
in  the  input.  Denote  the  input  value  M.  Then: 


if 

if 

if 


MS— I 
M5=  0 
MS=+I 


m  -  1 


M 

M  =  is 


MN* 

W7C 
MN1. 


MWC 


?  vv/c  -  r^ouif? 


where: 


C  Radial  pad  clearance,  inch 

D  Journal  diameter,  inch 

L  Bearing  length,  inch 

I  Mass  moment  of  inertia  of  pad^around 

longitudinal  axis,  lbs-in-sec 

2  2 
M  ■  I/R  equivalent  pad  mass,  lbs-sec  /in. 

N  Rotational  speed,  RPS 

R  Journal  radius,  inch 

uNDL  R  2 

S  »  ^ —  (^)  ,  overall  bearing  Sommerfeld  Number 

W  Total  bearing  reaction,  lbs. 

2 

U  Lubricant  viscosity,  lbs -sec/in 


In  general  it  is  convenient  to  use  the  last  form  of  M  (i.e.  for  MS  “  +1) 
since  it  is  independent  of  speed. 

Word  7  (NPR) .  If  NPR  =>  1,  the  output  includes  the  calculated  values  of 
the  pad  spring  and  damping  coefficients  in  addition  to  the  overall  bear¬ 
ing  coefficients.  If  NPR  =  0,  the  pad  coefficients  are  not  given,  only 
the  bearing  coefficients. 


List  of  Pad  Position  Angles 
4  •  (E15.7) 

This  input  list  defines  the  position  of  the  pads  by  means  of  the  angle  ^ 
measured  from  the  vertical  to  the  pivot  point  of  the  pad  (See  Fig. 22). 

Tj/  is  measured  In  degrees.  There  must  be  NP-values  (word  1,  card  3), 

4  values  per  card. 


List  of  Equivalent  Pad  Masses  if  MC  -  +1 
4(E15.7^ 


This  input  list  is  used  when  it  is  desired  to  study  the  effect  of  pad 
inertia  in  the  calculations  and  all  pads  have  identical  inertia.  This 
form  of  the  input  list  can  only  be  supplied  when  MC=+1  (word  5,  card  3). 
The  pad  inertia  is  given  in  form  of  the  dimensionless  equivalent  pad 
mass.M.  The  definition  of  M  depends  on  the  value  of  MS  (word  6,  card  3). 
The  input  list  must  contain  as  many  values  of  M  as  given  by  NM  (word  3, 
card  3),  4  values  per  card.  Thus  it  is  possible  to  obtain  the  dynamic 
bearing  coefficients  as  functions  of  the  pad  inertia. 


List  of  Pad  Forces  and  Derivatives 


The  spring  and  damping  coefficients  are  calculated  from  the  pad  forces 
and  their  derivatives.  These  quantities  can  be  computed  by  means  of 
PN0091,  Appendix  C,  and  data-reduced  by  program  PN0131,  Appendix  B. 

The  output  cards  from  PN0131  can  be  used  directly  to  make  up  the  present 
input  list.  The  input  list  consists  of  a  number  of  sets  of  complete  bearing 
data,  identified  by  a  bearing  eccentricity  ratio.  There  are  as  many  sets 
as  given  by  NEB  (word  2,  card  3).  Each  set  comprises  the  forces  and  their 
derivatives  for  each  pad.  To  illustrate: 


Card 

1  (3(E15 .7)  ) 

e. 

1st  Pad-< 

r 

Card 

2a  (4(E15.7)  ) 

e  « 

Card 

3a  (4(E15 .7)  ) 

afVde 

1st 

c 

^  Card 

4a  (4(E15 .7)  ) 

dftU 

0 

2nd  Pad  • 

j"  Card 

2b  (4(E15 . 7)  ) 

t  < 

p 

f* 


last  pad 


Card  4x  (4(E15.7)  ) 


2nd 


Card  1  ( 3 (El 5 . 7 )  ) 

Card  2  (4(E15.7)  ) 


1st  Pad 


where 


€c 


£ 

t 

i 

ft 

CO 


bearing  eccentricity  ratio,  see  Fig.  22 
bearing  attitude  angle,  degrees,  see  Fig.  22 
ratio  of  pivot  circle  clearance  to  pad  clearance 
pad  eccentricity  ratio,  see  Fig. 22 
pad  attitude  angle 

dimensionless  radial  force  component 
dimensionless  tangential  force  component 
angular  speed  of  journal,  rad/sec 


"Card  1"  above  is  actually  not  used  by  the  program  in  the  calculations. 

It  serves  only  the  purpose  of  identification.  Furthermore,  it  should  be 
noted  that  e  cannot  be  zero  since  the  calculations  include  f^/e  and  ^t/€‘ 


List  of  Equivalent  Pad  Masses  if  MC  ■  -1 
4(E15.7) 


This  input  is  used  when  the  pads  have  different  inertia.  This  form  of 
the  input  list  can  only  be  supplied  when  MC“-1  (word  5,  card  3).  The 
pad  inertia  is  given  in  the  form  of  the  dimensionless  equivalent  pad 
mass  M  as  defined  through  the  value  of  MS  (toord  6,  Card  3).  An  input 
list  must  be  given  for  each  bearing  eccentricity  ratio,  i.e.  following 
"Card  4x"  for  each  above.  The  list  is  made  up  of  NM  -  sets  (word  3, 
card  3)  and  each  set  contains  NP  values  of  M  (word  1,  card  3).  To 
illustrate: 


1st 


last  pad 


2nd  -T 
€o  1 


(■List 


1 

(  Card 

4x 

Kld(k 

pad 

1 

Mpad  2 

M  ,  - 
pad  3 

M  . 
pad 

5  “ 

M  pad  NP 

M  . 
pad 

1 

^pad  2 

~  -  -  - 

Mpad  NP 

Card 

1 

ft  c 

of  Pad  forces 
and  Derivatives 


M 


List  of  Equivalent 
V  Pad  Masses  for  MC  = 
NM  sets 


List  of  Pad  Forces 
and  Derivatives 
(continued) 
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OutPut  Data 

The  output  first  repeats  the  input  data  for  identification  and  checking 
purposes.  Then  follow  the  results  for  each  bearing  eccentricity  ratio. 

The  text  identifying  the  numerical  values  is  explained  below  except 
where  the  meaning  is  obvious: 

SOMMERF.NO.  -  When  part  of  pad  results:  S  =  (~)^  (F=total  force  on  pad) 

n  r  L 

When  part  of  bearing  results:  S=  (W*»total  force  on  bearing) 

F  -  T  =  f^ 

F  X  I  ■  -  1/Sp  ,  where  Sp  is  the  pad  Sommer feld  Number. 


K11,K12,K21,K22  =  K  gg  t  K^  (  i*e.  ^  spring  coefficients  for 

the  fixed  pad. 


WCll ,WCl2 ,WC21  ,WC22  =  ujCjj,  i.e.the  4  damping  coefficients 

for  the  fixed  pad. 

The  8  coefficients  are  dimensionless  in  the  form: 


-  f  ^rr  w  cff  *  ? 

where  F  is  the  total  force  on  the  pad.  Thus  Kll  is  cal 
culated  as: 


CRIT .MASS 


Kll  =  Sp[tfc«>s*<p  ~  ^  S inx(P  -($£  “  3e  -  t>] 


and  similarly  for  the  other  coefficients,  see  Eq.  (15). 
CFM 


crit 


(pDL(R}2f 


where  part  of  pad  results. 


CWM 


crit 


(pDL(R}2f 


where  part  of  bearing  results. 


KDll  ,KD12,KD21  ,KD22  =  KL,  Kl.  K'  Kl.  i.e.  the  spring  coefficients 

ff '  1 7 1  7F  /  77 

for  the  tilting  pad. 

WCD11  ,WCD12,WCD21  ,WCD22  =  (nCj^coC^ ,  cdC^  ,  ooC^  i.e.  the  danping 
coefficients  for  the  tilting  pad. 


The  coefficients  are  dimensionless  in  the  form: 


'  F 


Ff 


etc. 
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KXX,KXY,KYX,KYY  «  X  ,K  ,K  ,K  ,  i.e.  the  spring  coefficients  for 
xx  xy  yx  yy  6 

the  complete  bearing. 

WCXX ,  WCXY , WCYX , WCYY  =  03C  ,03C  ,coC  ,cuC  ,  i.e. the  damping  coefficients 

>  xx  xy  yx  yy  r  e 

for  the  complete  bearing. 

The  coefficients  are  dimensionless  in  the  form:  kix = 

OoC**  =  ^/CoCxx  etc.  where  W  is  the  total  force  on  the 
bearing.  Thus  for  ins  tance  is  calculated  internally 
in  the  program  directly  on  the  basis  of  the  input  values 
for  the  pad  forces  and  their  derivatives  and  the  final 
result  is  multiplied  by  the  bearing  Sommerfeld  Number  S. 


DIMENSIONLESS  PAD  MASS  -  is  the  pad  mass  as  defined  through  MS  (word  6, 
card  3,  input) . 

AMPLITUDE  *  \\[  i.e.  the  ratio  between  the  actual  pad  angular  amplitude 

and  the  angular  amplitude  of  a  massless  pad,  see  Eq.  (58). 
Hence,  the  ratio  may  be  thought  of  as  a  magnification  factor. 


PHASE  ANG  = 


CALC.  MASS 
MASS  1 
MASS  2 
MASS  3 


4r<j('y}0}~arj  (''Jp),  i.e.  the  phase  angle  in  degrees  between 
the  angular  motion  of  a  massless  pad  and  the  angular  motion 


of  the  actual  pad. 


1 

s  w/c 
i  mn£ 
s  w/c 

MN2 

w/c 

i  m£ 

s2  w/c 


MWC 


(^DL.R. 2)  J 
V 


Four  different  forms  of  the  dimenionlese 
pad  inertia.  For  details,  see  ex¬ 
planation  in  "Input  Data",  Card  3, 

Word  6. 


M  . 
CtU/M 


03 


crit 


/03, 


i.e. 


the  ratio  between  the  resonant 


pad  frequency  0Jcr^t  and  the  operating  angular  journal  speed 
for  the  given  Sommerfeld  Number. 


CRIT.  SP. RATIO 
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INPUT  FORM  FOR 

PN0078:  Spring  and  Damping  Coefficients  for  the  Tilting  Pad  Journal 
Bearing. _ _ 

IBM  1620  Computer.  FORTRAN  I 
Card  1;  Text,  Col.  2-52 


Card  2:  Text,  Col.  2-52 


Card  3:  7  •  (15) 


1.  NP.  Number  of  pads.  NP  6  12 

2.  NEB.  Number  of  bearing  eccentricity  ratios 

3.  NM.  Number  of  calculations  including  pad  inertia  per 

bearing  eccentricity  ratio.  NM  4  12. 


4.  NF. 

5.  MC. 

6.  MS 

6.  NPR 


If  NF  »  0, input  value  of  f  is  pad  Sommerfeld  number  and 

input  value  of  f  is  horizontal  pad  force  f*  .  If  NF  «■  1, 

f  =  f  and  f  =f  . 
r  r  t  t 

If  MC  =  -1:  each  pad  is  given  its  own  inertia 

If  MC  »  0:  all  pads  have  zero  inertia 

If  MC  =  +1:  all  pads  have  same  inertia  2 

If  MS  =  -1:  dimensionless  pad  inertia  =  s  ^C 

MN2 

If  MS  =  0:  dimensionless  pad  inertia  = 

_  _  1  MN2  MWC 

If  MS  =  +1:  dimensionless  pad  inertia  ■  ~2  ^j{ -  =  ^DL  R  2)2 

(q} 

If  NPR  =  0:  output  does  not  include  results  for  individual  pads 
If  NPR  =  1:  output  includes  results  for  individual  pads 


List  of  Pivot  Point  Angles:  4(E15.7) 

Give  the  position  angle,  degrees,  of  the  pivot  point  for  each  pad,  in  total 
NP  —  values,  4  values  per  card. 
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Llst  of  Dimensionless  Pad  Inertia.  MC  ■  +1:  4  •  (E15.7) 

Use  this  list  only  when  MC  *  +1.  Give  the  dimensionless  pad  inertia 
(same  for  all  pads)  according  to  value  of  MS,  in  total  NM  -  values, 

4  values  per  card: 


List  of  Bearing  and  Pad  Data:  4  (E15.7) 

Give  NEB-sets  of  cards.  If  MC  4  -1,  each  set  contains  (1+3  •  NP)  cards. 

If  MC  *»  -1,  each  set  contains  (1  +  3{fP  +  MN  •  cards: 

- £.  f.  C'/c 

- £  f  f,  ft  ) 

<jfr  if*  it  it 

-  -  -  -  Si  5e  3f 

_ */<)(£)  *74) . 

- £  f  f, 

it  it  &$r  it 

-  -  -  -  ie  si  iof  eif 

_ ‘S&ai  *&&  )  , 


—  ii  ■ 


—  H- 


—  u- 


U 

U 


—  ii 


—  n 


—  ii 


.4tl 


MP 


PJ 

^  inertia 
.VMCH 


C  H  C* 


2nd  Bearini 


fo  C'/C 

£  <P  <t 

&  &  &  dffr  ■  l4t. 

de  de  edf  tdf  pft« 

^/xj)  Jf'/a(i) 


£  f  (,  ft 

dfr  dft  &ft  M* 
3e  de  edf  edf 

dfM)  Mt/d(&) 


PU 

,  incrtiV 
If  MCH 


-  Mflu  |  ,  ,  ~  ~ 

-  NP 

-  R*di ,  "rug, 


MPaJ  np 
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DIMENSION  PVAI20) .CS(20) .SN<20) .CSCK20) »SNQ < 20 ) .CSN ( 20 ) .P < 20 > 
DIMENSION  B(8)»GXX(20) .GXY<20) .GYX<20) »GYY(20) .HXX < 20 ) »HXY < 20 ) 
DIMENSION  HYX ( 20 ) *HYY ( 20 ) .CMS(20) 

50  READ  200 
READ  201 

READ  202.NP.NEB.NM.NF.MC.MS.NPR 

PUNCH  203 

PUNCH  200 

PUNCH  201 

PUNCH  204 

PUNCH  205.NP .NEB.NM.NF.MC .MS.NPR 
PUNCH  206 
DO  51  1  =  1. NP. 4 

READ  20  7 »PVA ( I ) .PVA( 1  +  1  ) »PVA< 1  +  2) .PVA( 1+3) 

51  PUNCH  207.PVA(  I  )  .PVAU  +  l  )  .PVAU+2  )  *PVA<  1+3  ) 

IF  (MC)  54.54.52 

52  DO  53  1  =  1 .NM  .4 

53  READ  207.P( I  )  .P(  1  +  1  ).P<  1+2)  »P( 1+3) 

54  DO  55  1=1.8 

55  B ( I  )=0.0 

DO  100  1=1. NP 
C2=. 017453293*PVA( I  ) 

C 1 =COS { C  2 ) 

C2=S I N ( C2 ) 

CS(  I  )  =C  1 
SN ( I  )  =C2 
CSQ ( I )=C1*C1 
SNQ ( I )=C2*C2 

100  CSN ( I )=C2*C1 
DO  157  1=1. NEB 

READ  208.ECB.ATB.CRT 
PUNCH  210 

PUNCH  209.ECB.ATB.CRT 
BSN  =  0  •  0 
FH=0 • 0 

DO  120  J=1»NP 

READ  207.ECP.ATP.FR .FT. DFRE .OF TE.DFRA.DFT A. DFRS.DFTS 

C1=.017453293*ATP 

CSP=C0S(C1) 

SNP  =  S I N  <  Cl ) 

CSPQ=CSP*CSP 
SNPQ=SNP*SNP 
CS2=SNP*CSP 
IF(NF)  101.104.101 

101  C 1 =FR 
C2  =  FT 

FXI=-FR*CSP-FT*SNP 

ft=-fr*snp+ft*csp 

FR=SQRT(FR*FR+C2*C2) 

IF(FR) 103.102.103 

102  FR=1.0E90 
GO  TO  105 

103  FR=1.0/FR 
GO  TO  105 

104  FXI—1.0/FR 

105  BSN=BSN+FXI*CS( J)-FT*SN( J) 

FH=FH-FXI*SN( J)-FT*CS( J) 

PUNCH  211.  J 

PUNCH  212 

PUNCH  213. ECP.ATP.FR. FT. DFRE 

PUNCH  214  PN0078 
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IF  (NF)  106*1 07 ♦106 

106  PUNCH  215 

PUNCH  213. Cl *C2*FXI 

107  C1=(DFR£+DFTA)*CS2 
C2  ® ( OFT  E-DFRA )  *CS2 
C3=DFRS*CS2 
C4=DFT  S*CS2 
IF(ECP)  109*108*109 

108  C5=0.0 
C6=0.0 

GO  TO  110 

109  C5  =  FX I /ECP 
C6-FT /ECP 

110  AXX=DFRE*CSPQ-DFTA*SNPQ+C2-C6*SNP 
AXY=DFRA*CSPQ+DFTE*SNPQ+C1+C6*CSP 
AYX=-DFTE*CSPQ-DFRA*SNPQ+C1+C5*SNP 
AYY=-DFTA*CSPO+DFRE*SNPQ-C2-C5*CSP 
C5=2.0*C5 

C6=2.0*C6 

8XX=DFRS*CSPQ+C4-C5*SNP 
BXY=DFTS*SNPQ+C3+C5*CSP 
BYX=-DFTS*CSPQ+C3-C6*SNP 
BY Y=DFRS*SNPQ— C4+C6*CSP 
IF  ( AYY )  56*55.56 

55  00=0.0 
GO  TO  57 

56  C10=(AYY+BYY/AYY*BYY) /39.478418 

57  CMS ( J ) =C 10 

IF  ( NPR )  111.112.111 

111  C 1 =FR*AXX 
C2=FR*AXY 
C3=FR*AYX 
C4=FR*A YY 
C5=FR*BXX 
C6=FR*BXY 
C7=FR*BYX 
C8=FR*BYY 
C9=C10/FR 
PUNCH  216 

PUNCH  213. Cl .C2.C3.C4 
PUNCH  217 

PUNCH  213.  C5.C6.C7.C8»C9 

112  IF  (MC)  119,113.119 

113  C1-AYY*AYY+BYY*BYY 
IF  (Cl)  115,114,115 

114  PUNCH  218, J 
GO  TO  120 

115  C2=AXY#AYX-BXY*BYX 
C3=AXY*BYX+AYX*BXY 
C4=AXX-( AYY*C2+BYY*C3)/C1 
C5=BXX-(AYY*C3-BYY*C2)/C1 
C 1 *CSQ( J ) 

C2*SNQ( J ) 

C  3  =  CSN ( J ) 

B  ( 1 )  *B ( 1 ) +C4*C1 
B(2)«B(2)+C4*C3 
B  ( 3 ) *B ( 2 ) 

B(4)*B(4<+C4*C2 
B ( 5 ) *B( 5 )+C5*Cl 
B ( 6 ) =B ( 6 ) +C5*C3 


B(7)*B(6) 

B ( 8 ) =B ( 8 ) +C5*C2 
IF  (NPR)  117*118*117 

117  PUNCH  219.C4.C5 

118  GO  TO  120 

119  GXX ( J ) *AXX 
GXY  <  J ) » AXY 
GYX(J) =AYX 
GY Y ( J ) = AYY 
HXX ( J ) =BXX 
HXY ( J ) =BXY 
HYX ( J ) =BYX 
HY Y ( J ) =BYY 

120  CONTINUE 
6SN=1.0/BSN 

IF  ( MC )  123.121.123 

121  DO  60  J* 1 . 8 
60  B < J)=BSN*B( J ) 

PUNCH  220 

PUNCH  213»B(1)»B(2).B(3),B(4)»B(5) 
PUNCH  221 

PUNCH  213.B(6)  .8(7) .8(8 ) .BSN.FH 
PUNCH  228 
DO  158  J= 1 »NP 
C9=CMS( J)/BSN 
158  PUNCH  229,  J.C9 
DO  122  J»1.8 

122  B(J)=0.0 
GO  TO  157 

123  DO  156  K=1,NM 

IF  (MC)  124,126,126 
’24  DO  125  J- 1 ,NP .4 

125  READ  207. P(J)  ,P( J+l ) ,P( J+2) »P(J+3) 
GO  TO  127 

126  PM  =  P  (  K ) 

PUNCH  223.PM 

127  DO  151  J*1 »NP 

IF  (MC) 128,129,129 

128  PM=P ( J ) 

129  PZ=39.478418*PM 

IF  (MS)  130,131.132 

130  PC=PZ 
P  1  *PM 
P2=BSN*PM 
P3=PM/BSN 
GO  TO  133 

131  PC=PZ/BSN 
P 1 =PM/8SN 
P2  =  PM 
P3=P1/BSN 
GO  TO  133 

132  PC=PZ*BSN 
P 1 =PM*BSN 
P2=P1*BSN 
P  3  =  PM 

133  AXX=GXX ( J) 

AXY=GXY ( J ) 

AYX=GYX ( J ) 

AYY=GYY ( J) 

BXX=HXX ( J ) 

BXY=HXY ( J ) 
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BYX=HYX( J) 

B YY=HY Y ( J ) 

C3=AYY-PC 
C4=C3*C3+BYY*BYY 
IF  <C4)  137.134*137 

134  IF  ( AY Y )  136*135*136 

135  PUNCH  2 1 8 ♦  J 
GO  TO  151 

136  PUNCH  22 2  ♦  J 
GO  TO  151 

137  C1=C3/C4 
C2=BYY/C4 

C  3  =  1 • 0+C 1*PC 
C4=C2*PC 

ETP  =  SQRT  (C3*.C3+C4*C4) 

IF  ( C3 )  146*138*146 

138  IF  (BYY) 142*139.142 

139  IF  ( AYX+BYX )  141*140*141 

140  PUNCH  218. J 
GO  TO  151 

141  ETP=0.0 
PHA=180.0 
GO  TO  150 

142  IF  (C4)  145*143*144 

143  PHA=0.0 
GO  TO  150 

144  PHA=90 • 0 
GO  TO  150 

145  PHA=-90 • 0 
GO  TO  150 

146  C4=C4/C3 
PHA=57.295780*ATAN(C4) 

IF  ( C 3  >  147*150*150 

147  IF  <C4)  148.150.149 

148  PHA=PHA+180.0 
GO  TO  150 

149  PHA=PHA— 180*0 

150  C5=-PC*C3 
C6=PC*C4 

C7=C1*AXY+C2*BXY 

C8=C1*BXY-C2*AXY 

GR=AXX-C7*AYX+C8*BYX 

GE=BXX-C7*BYX-C8*AYX 

DR=~PC*C7 

DE=-PC*C8 

C3=_PC«(C1*AYX+C2*BYX ) 
C4=-PC*(C1*BYX-C2*AYX) 
C8=CMS( J) /BSN 
IF  ( P  3 )  161*160.161 

160  C9= 1 *0E+09 
GO  TO  162 

161  C9=C8/P3 

162  PUNCH  211*  J 
PUNCH  224 

PUNCH  2 13 *GR *DR *C3 »C5 »GE 
PUNCH  225 

PUNCH  213*DE»C4»C6»ETP*PHA 
PUNCH  226 

PUNCH  2 13 »PC *P1 »P2 *P3 

PUNCH  227 

PUNCH  213*  C8.C9 
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C8=CSN( J) 

Cl= (DR+C3)*C8 
C2=  ( DE+C4 ) *C8 
C  7= ( GR-C5 ) *C  8 
C8  = ( GE-C6 ) *C8 
AXX=CSO ( J) 

8XX  =  SNQ ( J ) 

B ( 1  )  -B< 1 )+GR*AXX+C5*BXX-Cl 
B(2  )  =B ( 2  )+DR*AXX-C3*BXX+C7 
B  (  3  )  =B ( 3)+C3*AXX-DR*BXX+C7 
B(4)=B(4)+C5*AXX+GR*BXX+C1 
B<5  )=B( 5>+GE*AXX+C6*BXX-C2 
B(6)=B(6)+DE*AXX-C4*BXX+C8 
B(7)*B< 7)+C4*AXX-DE*BXX+C8 
B(8)=B( 8)+C6»AXX+GE*BXX+C2 

151  CONTINUE 

IF  ( BSN  )  152*154.152 

152  DO  153  J= 1 » 8 

153  B ( J ) sBSN*B ( J  ) 

154  PUNCH  220 

PUNCH  213»B(1).B(2)»B(3)»B(4)»B(5) 

PUNCH  221 

PUNCH  213. B<6) *B(7)  •  B  <  8 ) .BSN.FH 
DO  155  J= 1 » 8 

155  B(J)=0.0 

156  CONTINUE 

157  CONTINUE 
STOP 

200  FORMAT ( 49H0  ,3H  ) 

201  FORMAT ( 49H0  ,3H  ) 

202  FORMAT!  15.15. 15. 15, 15. 15. 15) 

203  FORMAT (49H1PN0078  SPR-DAMP .COEFF . FOR  TILT. PAD  JRNL .BRG. 1-18 , 3H-63 ) 

204  FORMAT ( 49H0  N.PDS  N.ECC  N.MSS  FORCE  M.INP  M.TPE  PRT ) 

205  FORMAT (1X15,2X15,2X15,2X15.2X15.2X15,2X15) 

206  FORMAT ( 19HOPIVOT  POINT  ANGLES) 

20  7  FORMAT  <E15.7,E15.7,E15.7,E15.7) 

208  FORMAT (E15.7.E15.7.E15.7) 

209  FORMAT ( 9H  BRG . ECC= » E 14 . 7 , 1 2H  ATT.ANG= ,E14.7 ,9H  CP/C=.E14.7) 

210  FORMAT ( 2HQ  /20H0BEAR I NG  CALCULATION) 

211  FORMAT ( 11HOPAD  NUMBER, 13) 

212  FORMAT (4X8HECC. RAT. 7X7HATT.ANG5X10HSOMMERF.NO8X3HF-T9X6HDFR/DE ) 

213  FORMAT ( 2X  E12.5.2X  E12.5.2X  E12.5.2X  E12.5.2X  E12.5) 

214  FORMAT ( 5X6HDFT /DE8X7HDFR/EDA7X7HDFT / EDA6X8HDFR/DE/ W6X8HDFT / DE / W ) 

215  FORMAT ( 7X2HFR 12X2HFT 1 2X3HFX I  ) 

216  FORMAT < 7X3HK1 1 1 1 X3HK 1 2 1 1X3HK2 11 IX 3HK22) 

217  FORMAT ( 6X4HWC 1 1 10X4HWC 1 2 1 0X4HWC2 1 10X4HWC228X9HCR I T . MASS ) 

218  FORMAT ( 18H0M0TI0N  OF  PAD  NO.,l3»l4H  INDETERMINATE) 

219  FORMAT ( 6H  KD1 1 =E 12 . 5 , 6X6HWCD1 1 = E 1 2 . 5 ) 

220  FORMAT  (  2H0  /1H  6X3HXXX 1 1 X 3HK.XY  1 IX 3HX YX1 1  X3HK.Y Y 10X4HWCXX  ) 

221  FORMAT (6X4HWCXY 10X4HWCYX1 0X4HWCYY 6X1 1HSOMMERF.NO. 4X1 OHHOR 12. FRC.  ) 

222  FORMAT ( 33HORESONANCE  WITHOUT  DAMP  I NG  »  PAD  NO.I3) 

223  FORMAT ( 24H0D I  MENS  I ONLESS  PAD  MASS=»E12.5) 

224  FORMAT (6X4HKD11 10X4HKD1 21 0X4HKD21 10X4HKD221 0X5HWCD11  ) 

225  FORMAT (6X5HWCD129X 5HWCD219X 5HWCD227X9HAMPL I TUDE5X9HPHASE  ANG) 

226  FORMAT ( 4X9HCALC.MASS6X6HMASS  18X6HMASS  28X6HMASS  3) 

227  FORMAT (4X9HCR IT. MASS2X13HCRI T.SP.RAT 10) 

228  FORMAT (9X7HPAD  NO. 2X9HCR I T .MASS ) 

229  FORMAT (8X15, 3XE12. 5) 


END 
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COMPUTER  PROGRAM  PN0131:  INTERPOLATION  OF 
PARTIAL  ARC  BEARING  FORCE  DERIVATIVES  FOR 
USE  IN  TILTING  PAD  BEARING  CALCULATIONS 
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APPENDIX  B 

Computer  Program  PN0131:  Interpolation  of  Partial  Arc  Bearing 
Force  Derivatives  for  Use  in  Tilting  Pad  Bearing  Calculations 

In  order  to  calculate  the  spring  and  damping  coefficients  for  the 
tilting  pad  journal  bearing  it  is  necessary  to  know  the  force  deri¬ 
vatives  of  the  pads  making  up  the  complete  bearing.  Since  these 
derivatives  are  calculated  numerically  in  a  computer  by  finite 
difference  methods  they  are  usually  only  available  for  a  limited 
number  of  eccentricity  ratios.  When  the  pad  is  incorporated  in 
the  tilting  pad  bearing  its  eccentricity  ratio  is  determined  by 
the  bearing  eccentricity  ratio  (see  Eq.  (B6))  and  it,  therefore? 
becomes  necessary  to  know  the  force  derivatives  at  other  eccentri¬ 
city  ratios  than  the  ones  for  which  calculations  have  been  per¬ 
formed.  It  is  rather  costly  and  impractical  to  obtain  this  data 
by  performing  additional  calculations  with  the  finite  difference 
computer  program  and  instead  interpolation  of  the  existing  data 
can  be  employed.  It  has  been  found  that  to  do  this  manually  (by 
curve  plotting)  is  time  consuming.  Hence,  a  computer  program  has 
been  written  for  the  purpose  of  performing  the  interpolation.  The 
computer  program  accepts  as  input  a  list  of  the  force  derivatives 
as  computed  from  the  finite  difference  program  (i.e.PN0091,  Appendix  C). 
In  addition  the  number  of  pads  and  their  position  must  be  given.  Then 
the  computer  program  calculates  the  force  derivatives  for  each  pad  for 
a  specified  number  of  bearing  eccentricity  ratios.  The  output  can  be 
used  directly  as  input  for  the  tilting  pad  bearing  program  (i.e.  PN0078, 
Appendix  A) . 

Analysis 

The  interpolation  is  based  on  quadratic  curve  fitting.  Let  the  function 
to  be  interpolated  be  denoted  0^  which  is  given  for  a  number  of  discreet 
eccentricity  ratio  values  by  the  input.  Then: 

Ai (e-e if +'Bj  (e-et)  +  ^  5-,  ^ 


(Bl) 
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where: 

a  =  ifru.-ft)- V«.-o.-) 

(B2)  A,+A2 

(B3)  B  =  ^4,-., -$.■)•*■  I? 


^2  r 


■l+l 


Thus,  for  each  £-  interval,  except  the  first  and  the  last,  ^  is 
determined  by  two  equations.  The  average  value  is  used: 

<B4)  +  A  £•£;.,  f+B  J  +  ft.,  ] 


There  are  nine  functions  to  be  interpolated.  Since  they  all  (except 
the  attitude  angle)  become  infinite  when  £=!  it  is  found  necessary 
to  normalize  the  functions  before  interpolation  takes  place.  The 
normalization  is  chosen  as  follows: 


u± 


Wo  norma  liZAtich 


(B5) 


€l£  £  -  6h 


r,  is  d'M  h  0-e1)* 

Wor*h*lr2At.(»n  4<ecV:  IS  S  fl-ef'* 

&  if*  •  j.  j  j  l 

£d| ;  idf  w  h  l-e* 


Here  £,  and  £h  are  the  first  and  last  eccentricity  ratio,  respectively, 
in  the  input  list,  flu  is  specified  by  the  input  as  described  later. 


Let  the  tilting  pad  bearing  have  In  pads  whose  pivot  points  are  located 
the  angle  from  the  vertical,  measured  in  the  direction  of  rotation. 
Then  the  eccentricity  ratio  of  any  pad  is  determined  by: 


(B6) 


tcosf*  l-  i' -  £tcos(y-cpc ) 


where  £0  is  the  bearing  eccentricity  ratio,  fo  is  the  bearing 
attitude  angle,  (J-f  j  is  the  preload,  £  is  the  pad  eccentricity  ratio  and 
< p  is  the  pad  attitude  angle.  The  input  specifies  flnd  y 

and  (ft  is  given  as  a  function  of  £  in  the  earlier  mentioned  input  list. 


Thus,  ECOSf  can  be  calculated  from  the  Input  list  and  by  Inter¬ 
polation  Eq.  (B6)  can  be  solved  to  find  £  for  each  pad.  With  £ 
determined  the  8  force  parameters  can  be  calculated  by  Interpolation 
using  Eq.(B4)  and  (B5). 

Input  Data 

The  program  is  written  for  the  IBM  1620  computer,  40K  memory  storage, 
with  Input  and  output  punched  on  cards. 

Card  1 

Descriptive  text,  Column  2  to  52 

Card  2 
5  •  (1X14) 

Word  1  (NEP)  gives  the  number  of  pad  eccentricity  ratios  for  which  the 
force  derivatives  are  provided  in  the  pad  data  list.  NEP  ^  15. 

Word  2  (LM)  gives  the  value  of  L  as  used  in  Eq.(B5).  If  the  pad 
eccentricity  ratios  in  the  pad  data  list  are  £h  then  the 

pad  data  are  divided  by  a  scale  factor  internally  in  the  program  for 
£  ^  £l  •  Since  the  scale  factor  is  zero  for  £.-0  it  has  been  found 
necessary  to  set  the  scale  factor  equal  to  1  for  £*£L,  In  the 
determination  of  L  the  most  critical  parameter  is  fsf  which  is  usually 
negative  for  .5,  Hence,  it  is  suggested  to  give  L  such  a  value 

that  £W|sS',5‘  +o  .6  • 

Word  3  (NBE)  gives  the  number  of  tilting  pad  bearing  eccentricity  ratios. 
There  is  no  limit  to  NBE. 

Word  4  (NA)  gives  the  number  of  pads.  NA  ^  20. 

Word  5  (INP)  If  INP=0,  more  input  data  follows  the  present  set  of  input. 


If  INP  4  0,  the  present  set  of  input  is  the  last  set. 
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Pad  Data  List 
5  ■  (1XE13.6) 

The  pad  data  is  usually  obtained  from  a  finite  difference  solution  of 
Reynolds  equation,  see  Appendix  C.  The  pad  forces,  force  derivatives 
and  attitude  angle  are  calculated  for  a  number  of  pad  eccentricity  ratios 

£i  , fj - £h  and  supplied  as  input  to  the  present  computer  program.  For 

each  eccentricity  ratio  two  cards  are  given: 


First  card:  £ 

Second  card:  4+ 


Ms 

£d<f 


where: 

£  =  pad  eccentricity  ratio 

<p  =  pad  attitude  angle 
,  Fr 

V  “  yuN/DL^/fc)1  j  dimensionless  radial  force 
,  Ft 

=  jn  NDL(R/c)*  dimensionless  tangential  force 
£  =  radial  journal  center  velocity,  sec  * 


Here,  (f)j  fr  and  ft  are  calculated  for  the  given  £  at  the  steady  state 
position  such  that  the  total  force  passes  through  the  pivot  point. 


Hence,  the  pad  data  list  consists  of  2  * (NEP)  cards  (see  word  1,  card  2). 

The  list  is  given  in  sequence  such  that  £I<-£1£ - <  ^Ner 

List  of  Pivot  Point  Angles 

5  •,„(UE.U.i6) 

This  input  list  gives  the  values  of  the  angle  measured  from  the  vertical 
to  the  pivot  point  of  each  pad  in  the  direction  of  rotation,  see  Fig.  22. 

The  angle  is  measured  in  degrees.  There  are  (NA) -values  in  total  (word  4, 
card  2),  5  values  per  card. 


List  of  Bearing  Eccentricity  Ratios 
3  •  (E15.7) 


4S 


The  position  of  the  journal  center  with  respect  to  the  bearing  center 
is  defined  through  the  eccentricity  ratio  and  the  attitude  angle  <p0 
see  Fig.  22.  £0  is  calculated  with  respect  to  the  pad  clearance 

(pad  radius  minus  journal  radius)  and  is  in  degrees,  measured 

from  the  vertical  to  the  line  of  centers  in  the  direction  of  rotation. 
In  the  general  case  (off-center  pivot  position,  unsymmetric  arrange¬ 
ment  of  the  pads)  a  considerable  effort  is  involved  in  determining  (p0 
but  if  the  pads  are  centrally  pivoted  and  located  symmetrically  around 
the  vertical  then  ^,-0. 

There  must  be  an  input  card  for  each  bearing  eccentricity  ratio,  in 
total  (NBE)  cards.  Each  card  contains  three  words: 

Word  1  gives  the  bearing  eccentricity  ratio  £a 

Word  2  gives  the  bearing  attitude  angle  (f>0  in  degrees 

Word  3  gives  the  ratio  C/c  where  C‘  is  the  radius  of  the  pivot 

point  circle  minus  the  journal  radius  and  C  is  the  pad  radius  minus 

the  journal  radius. 

Output  Data 

The  output  gives  first  the  values  of  the  input  data  for  identification 
and  checking  purposes.  Then  follow  6  lines  of  text  describing  the 
format  of  the  output  data  and  thereafter  the  numerical  values  of  the 
output  data  arranged  as  follows  (FORMAT  4(E15.7): 


Card 

1,  2  and  3 

Card 

4 

-  c'fc. 

Card 

la 

- 

Card 

2a 

-  <fA  fr,  ft 

■  1st 

pad 

Card 

Card 

3a 

4a 

-*%$!*/*  f 

Card 

Card 

lb 

2b 

i  i 

rn^ 

n» 

>  2nd 

pad 

I 


I 

1 

I 


I 


Note  that  two  values  of  the  pad  attitude  angle  (f  are  given.  The 
first  value  (i.e.  on  card  la,  lb  etc)  is  obtained  by  interpolation 
of  the  pad  data  list  whereas  the  second  value  (i.e.  on  card  2a, 2b  etc) 
is  calculated  from  (j P-tan'(f*/fr) ,  It  is  the  latter  value  that  should 
be  actually  used  and  by  comparison  with  the  first  (jP-value  some  in¬ 
dications  are  obtained  of  the  accuracy  of  the  interpolation. 

The  output  cards  can  be  used  directly  as  input  to  the  tilting  pad  bearing 
computer  program  (PN0078,  Appendix  A)  by  removing  the  3  blank  cards  in 
front  of  each  bearing  eccentricity  ratio  calculation  (card  1,2  and  3  above) 
and  by  removing  the  first  card  of  each  pad  data  output  set  (card  la, lb  etc. 
above) . 

Returning  to  Eq.  (B6)  it  can  be  seen  that  geometrical  considerations  limit 
the  range  of  the  bearing  eccentricity  ratio  £  .  The  program  takes  care 

of  this  automatically  as  follows.  For  a  given  and  the  right 

hand  side  of  Eq.  (B6)  may  be  calculated.  Denote  the  result  . 

Three  cases  may  arise: 

a)  If  e,CoS(f^U(cS(f)(^l  y  Eq.  (B6)  is  solved  by  interpolation.  Here  £, 
is  the  first  eccentricity  ratio  value  in  the  pad  data  list. 

b)  If  (tcoS<f)CllU  4.  €,(0S(f,  the  computer  writes:  "NO  LOAD,  PVT.  ANG=X.XX, 

E  *  COS(A)=X.XX"  indicating  that  the  particular  pad,  identified  by  its 
pivot  angle  has  a  value  less  than  the  lowest  value  in  the 

pad  data  list.  Strictly  speaking  there  may  be  solutions  in  this  case 
but  assuming  that  € ,  in  general  is  almost  zero  it  is  safest  not  to  allow 
extrapolation. 

c)  If  (t(oS(p)f£l  the  computer  writes:  "E=l,  PVT.  ANG=X.XX,  E  *  COS(A)=X.XX" 
indicating  that  for  the  particular  pad  the  journal  is  either  touching 

or  beyond  the  pad  surface.  The  calculations  proceed  with  the  next  pad. 


-Mr- 


INPUT  FORM  FOR 

PN013I:  Interpolation  of  Partial  Arc  Bearing  Force  Derivatives  for 
_ Use  in  Tilting  Pad  Bearing  Calculations _ 

IBM  1620  Computer.  FORTRAN  I 

Card  1  (Text  Col. 2-52) 


Card  2:  5(1X14) 


1. 

NEP: 

Number  of  pad  eccentricity  ratios  in  pad 
list.  NEP  4  15 

data 

2. 

LM: 

=  pad  data  not  normalized! 

=  pad  data  is  normalized  J 

Eq.  (B5) 

3.  NBE: 

Number  of  bearing  eccentricity  ratios 

4. 

NA: 

Number  of  pads.NA  4  20 

5.  INP: 

If  INP=0:  more  input  follows.  If  INP^O, 
set  of  input. 

last 

Pad  Data  List:  5(1XE13.6) 


Give  2  cards  per  pad  eccentricity  ratio,  in  total  2  •  NEP  cards: 


e, 


—  II  - 


—  ll  — 


—  H  — 


— H  — 


— li  — 


-U- 


-  li  — 
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Plvot  Point  Angles:  5(1XE13.6) 

Give  the  position  angle  (degrees) of  the  pivot  point  for  each  pad, 
in  total  NA-values ,  5  values  per  card: 


Bearing  Eccentricity  Ratio:  3(E15.7) 

Give  one  card  per  bearing  eccentricity  ratio,  in  total  NBE  cards: 

- -  ec  ,  <p0  t  c'/c 


DIMENSION  P( 11*15) »PN  I 11*15) »A  ( 1 1  ♦  1 5  ) »B( 11.15). PA (20)»QI 11 )  46 

49  READ  101 

READ  102*NEP*LM»NBE»NA*INP 
PUNCH  100 
PUNCH  101 
PUNCH  104 

PUNCH  103.NEP.LM.NBE.NA. INP 
PUNCH  107 
DO  47  1=1.20 

47  PA (  I  ) =0 .0 
DO  48  1=1*11 
0 ( I ) =0  •  0 
DO  48  J  =  1  *  1 5 
A ( I ,J>*0.0 
B ( I »J)=0.0 
P( I  .J)=0.0 

48  PN ( I.J)=0.0 
KM=LM+1 
NM=NEP— 1 
DO  55  J= 1  * NEP 
PUNCH  108. J 
DO  50  1=1. 6. 5 

READ  109, P( I »J) »P( 1  +  1 »J) »P( 1+2. J) »P( I +3  *  J ) .P( I +4  ,  J ) 

50  PUNCH  109. P( I .J) .P( I+l.J) »P( l +2 ♦ J ) »P ( l  +  3  » J ) »P( I +4 » J ) 

C1=P( 1 .J) 

C2=P(2.J) 

C2=.017453293*C2 
C2=COSF (C2 I 
P(11»J)=C1*C2 
IF(J-LM)  55,51,51 

51  C2=1.0-C1*C1 
C3=C1/C2 
C4=SQRTF(C2) 

C4=C3/C4 
C2=C3*C3 
DO  52  1=3.5 
PN ( I.J)=P( I,J)/C2 

52  PN ( I +3  »  J ) =P ( I+3.J)/C4 
PN(9,J)=P(9,J)/C3 
PN(10.J)=P(10,J)/C3 
PN(2»J)=P(2, J) 

55  CONTINUE 
PUNCH  110 
DO  56  I =1 »NA  ♦  5 

READ  109 .PA I  I  ) ,PA( 1  +  1 ) ,PA< 1  +  2 ) »PA( 1+3) »PA( 1+4) 

56  PUNCH  109.PA ( I ) »PA( 1+1) ,PA( 1+2) »PA( 1+3) ,PA( 1+4) 

PUNCH  113 
PUNCH  114 
PUNCH  115 
PUNCH  116 
PUNCH  117 
PUNCH  118 
DO  65  1=2.10 
Y2  =  P(  I  .1 ) 

Y  3=P ( 1.2) 

X2*P( 1  ♦  1 ) 

X3=P( 1,2) 

DO  65  J=2.NM 
Y1  =  Y2 
Y2  =  Y3 
X  1  =X2 


PN0131 


X2  =  X3 

X3=P( 1 tJ+1  ) 

1  F  (  J-K.M )  60*61*62 

60  Y3®P< I *J+1 ) 

GO  TO  64 

61  Y 1  =PN ( 1 *LM) 

Y  2  =PN ( I »KM ) 

62  Y 3  =  PN ( I ♦ J+l ) 

64  DT1=X2-X1. 

DT2=X3-X2 
C1=Y2-Y1 
C2=Y3-Y2 
C  3  =  X  3-X 1 
C4=DT1/0T2 

A  < I .J)=<C2*C4-C1 ) / <  DT1*C3  » 

B ( I »J)®(C2*C4+C1/C4)/C3 

65  CONTINUE 
Y2=P( 1*1) 

Y3*P< 1.2) 

X2  =  P( 11 . 1 ) 

X3*P( 11.2) 

DO  66  J  =  2  »NM 
Y1  =  Y2 

Y2  =  Y3 

Y  3  =  P I  1 »J  +  1 ) 

X 1  =  X2 

X2  =  X3 

X3=P( ll.J+1 ) 

DTI =X2-X 1 

DT2=X3-X2 

C 1 = Y2-Y 1 

C2=Y3-Y2 

C3=X3-X1 

C4=DT1/DT2 

A ( 1 ,J) - (C4*C2~C1 ) / <  DT1*C3 ) 

66  B(1»J)=(C4*C2+C1/C4)/C3 
DO  239  K  =  1 » NBE 

READ  105  *EB *  ATB  »CPC 
PUNCH  112 

PUNCH  105.EB.ATB.CPC 
ATR=.017453293*AT8 
DO  238  L  =  1 »  NA 
PAN=PA ( L ) 

PAR®.017453293*PAN 
ECS® COS F ( PAR-ATR ) 

ECS® 1 • 0-CPC-EB*ECS 
I F  C  P ( X 1  ♦  1 ) -ECS  )  202*205*201 
201  PUNCH  1 19 » L  .PAN . ECS 
GO  TO  238 

102  I F ( 1 • Q-ECS )  203.203.204 

203  PUNCH  120*L.PAN*ECS 
GO  TO  238 

204  I F ( ECS-P ( 1 1 »NEP ) )  206.200.200 
200  KC=NEP 

GO  TO  209 

205  KC=2 

GO  TO  209 

206  DO  208  J®2  *NEP 
C1=P( ll.J) 

IF(ECS-Cl)  207,207.208 

207  KC  =  J 
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go  TO  209 

208  CONTINUE 

209  AR  =  A< l.KC) 

BR  =  B( 1 »KC ) 

CR=Pl l.KC) 

XR  =  P  (  1 1 » KC ) 

KL=KC-1 

AL  =  A  ( 1 ,KL ) 

BL  =  B<  1  » KL  ) 

CL  =  P( l.KL) 

XL  =  P< 11. KL) 

C  1  =ECS-XR 

C1=CR+C1*(BR+C1*AR) 

C2=ECS-XL 

C2=CL+C2*(bL+C2*AL) 

I F ( KC-2  >  210.210.211 

210  C2  =  C1 

GO  TO  213 

211  IF(NEP-KC)  212.212.213 

212  C 1 =C2 

213  EP=(Cl+C2)/2.0 
C2*1.0-EP*EP 
C3*EP/C2 
C4=SQRTF (C2 ) 

C4=C3/C4 

C2=C3*C3 

DO  234  1=2.10 
AL=A ( I . KL ) 

BL=B ( I .KL> 

XL=P( 1 .KL) 

AR=A ( I »KC> 

BR  =  B ( I » KC ) 

XR=P( l.KC) 

IF  (1-2)  220,220.221 
220  0  =  1.0 

GO  TO  226 

221  IF  ( 1-5  )  222  .222,223 

222  C 1 =C2 

GO  TO  226 

223  IF  (  1-8 )  224.224,225 

224  C 1 =C4 

GO  TO  226 

225  C 1  =C 3 

226  QR=EP-XR 
QR=QR*(BR+QR*AR) 
QL=EP-XL 
QL=QL*(BL+QL*AL) 

IF  (KC-LM)  227,227,230 

227  QR=QR+P (l.KC) 

228  QL=OL+P< I »KL ) 

IF  (KC-2)  229,229.233 

229  QL=QR 

GO  TO  233 

230  0R=C1* ( QR+PN (l.KC)) 

IF  (KC-KM)  228,228,231 

231  QL  =  C 1* ( QL+PN ( I  ,KL )  ) 

IF  ( KC-NEP )  233,232,232 

232  QR=QL 

233  Q ( I ) = ( OR+QL ) /2.0 

234  CONTINUE 
C2  =  0  (  3  ) 
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IF  (C 2)  236.235.236 

235  C  1  =0 ( 2 ) 

GO  TO  237 

236  C1=ATANF(Q(6)/C2) 

C1*57.295780*C1 

237  PUNCH  111,  PAN  »ECS  »Q<  2  ) 

PUNCH  106.  EP.C1.C2.Q<6) 

PUNCH  106.  Q(4) »Q< 7) ,Q<9)  ,Q< 10) 

PUNCH  106.  0 ( 5 ) »Q ( 8  > 

238  CONTINUE 

239  CONTINUE 

IF  (INP)  2 40,49,240 

240  STOP 

100  FORMAT (49H1PN0131-INTERPOLAT ION  OF  PART. ARC  DATA  J.LUND  8-9.3H-63) 

101  FORMAT ( 49H0  » 3H  ) 

102  FORMAT < 1X14. IX 14 .1X14, IX 14, IX 14) 

103  FORMAT ( 2X14 .6X14, 6X1 4.6X14, 6X1 4) 

104  FORMAT ( 9H0N0. PAD  E4X5HL I M I T3X8HNO • BRG . E3X7HNO. ANG. 5X5H I NPUT ) 

105  FORMAT (E15.7.E15.7.E15.7) 

106  FORMAT ( El 5 . 7 , El  5 . 7 , El  5 . 7 . E 1 5 . 7 ) 

107  FORMAT ( 49H0E  ATT  FR  DFRDE  DFRDEW  FT  DFTDE  DFTDEW  DFREDA  DFT.3HEDA) 

108  FORMAT ( 3H0E ( ♦ I  2 . 1H ) ) 

109  FORMAT  <1XE13.6,1XE13.6,1XE13.6,1XE13.6,1XE13.6) 

110  FORMAT ( 13H0PIVOT  ANGLES) 

111  FORMAT (1H0.E14.7.E15.7.E15.7) 

112  FORMAT ( //1H0 ) 

113  FORMAT ( 14H00UTPUT  FORMAT) 

114  FORMAT ( 8H  EPS.BRG5X7HAtT.BRG5X4HCP/C ) 

115  FORMAT ( 38H  PIVOT  ANG  E*COS(A)  ATT .PAD (  I NTP )  ) 

116  FORMAT ( 8H  EPS.PAD5X7HATT.PAD5X2HFR8X2HFT ) 

117  FORMAT ( 7H  DFR/DE6X6HDFT /DE6X7HDFR/EDA3X7HDFT /EDA ) 

118  FORMAT ( 23H  DFR/D( E/W)  DFT/D ( E/W)  ) 

119  FORMAT ( 4H0PAD . I  2  » 1 7H  NO  LOAD »PVT. ANG* »E1 3.6 » 10H , E*COS C A ) = ,E 1 3.6 ) 

120  FORMAT (4H0PAD.I2.13H  E  =  1 , PVT . ANG*  .El  3 .6 » 10H »E*COS ( A ) = ,E 1 3.6 ) 

END 
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A8?IMCT 

This  report  describes  the  computer  progress  PN0091: 
"Hydrodynamic  Incompressible  Journal  Bearing."  The 
program  exists  in  2  versions,  one  for  the  IBM  704 
computer  and  one  for  the  IBM  1620  computer.  The  re¬ 
port  contains  a  description  of  the  analysis  and  the 
instructions  for  using  the  programs.  The  program 
solves  the  dynamical  incompressible  Reynolds  equation 
in  dimensionless  form  by  finite  difference  equations 
and  gives  results  for  the  pressure  distribution,  the 
load  and  the  attitude  angle.  The  rupture  of  the  oil 
film  is  included.  The  bearing  may  be  full  360°  or  a 
partial  arc. 


INTRODUCTION 


The  dynamical  Reynolds  equation  is  a  partial  differential  equation  re¬ 
lating  the  pressure  in  the  fluid  film  to  the  journal  rotational  speed 
and  the  squeeze  film  velocity.  The  film  is  assumed  to  be  isothermal 
so  that  the  viscosity  is  constant.  The  equation  is  used  in  dimension¬ 
less  form  and  is  given  as  Eq.  (4)  page  4. 

The  bearing  is  taken  to  be  operating  in  air.  Therefore,  no  oil  can 
enter  the  bearing  from  the  ends  and  the  film  pressure  cannot  be  lower 
than  ambient.  For  this  reason  the  film  contracts  in  the  diverging  part 
of  the  bearing  with  a  uniform  pressure  of  .  This  is  denoted  film 
rupture  and  is  taken  into  account. 

Reynolds  equation  is  transformed  into  finite  difference  form  and  solved 
by  iteration.  The  convergence  of  the  iteration  process  is  established 
in  two  ways,  one  in  which  the  relative  difference  between  two  iterations 
is  checked  against  a  preassigned  error  and  one  which  compares  the  result 
of  each  iteration  against  an  extrapolated  value,  i.e.s  a  form  of  absolute 
convergence. 

The  program  provides  for  the  calculation  of  force  derivatives  with  respect 
to  eccentricity  ratio  and  the  two  squeeze  film  velocity  components.  Such 
derivatives  are  needed  in  computing  spring  and  damping  coefficients  for  the 
bearing. 

If  desired  the  program  will  iterate  on  the  attitude  angle  in  order  to  de¬ 
termine  the  attitude  angle  corresponding  to  no  horizontal  force  component. 
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analysis 

NQMBNCLAHJEE 

C  Diametrical  clearance,  inch 

D  Journal  diameter,  inch 

F  Force,  lbs. 

fr,ft  Radial  and  tangential  force  components,  dimensionless,  see  Figure  l. 
fy,fh  Vertical  and  horizontal  force  components,  dimensionless,  see  Figure  1 
h  Dimensionless  film  thickness 

i,j  Finite  difference  coordinates,  axial  and  circumferential,  see  Figure  2. 
k  Iteration  number 

L  Bearing  length,  inch 

m  Number  of  circumferential  subdivisions  see  Figure  2 

n  Number  of  axial  subdivisions,  see  Figure  2 

N  Journal  speed,  RPS 

P  Dimensionless  pressure  (above  ambient) 

R  Journal  radius,  inch 

5  Sommerfeld  number 

t  Time,  seconds 

U^,U2  Journal  and  bearing  surface  speed,  inch/sec. 

x,z  Circumferential  and  axial  coordinates,  dimensionless 

yfc  Sum  of  all  pressures  after  the  k'th  iteration,  see  Eq.  (10) 

y„  Extrapolated  pressure  sum  after  the  k'th  iteration,  see  Eq.  (11) 

Ct  Attitude  angle,  degrees 

a  Tangential  speed  of  journal  center,  rad/sec. 

6^  Absolute  convergence  limit,  see  Eq.  (12) 

6  Relative  convergence  limit,  see  Eq.  (9) 

&A  Error  in  absolute  convergence,  see  Eq.  (12) 

Error  in  relative  convergence,  see  Eq.  (9) 

g  Eccentricity  ratio 

€  Radial  speed  of  journal  center,  sec 

0  Circumferential  angular  coordinate,  see  Figure  1. 

2 

U  Viscosity,  lbs-sec/in 

O)  Angular  speed  of  journal,  rad/sec. 
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For  an  incompressible  lubricant  Reynolds  equation  is: 

<»  Mi£$f]  +  &[l£$l]  =iMK(U,+U,)]+  If 


To  make  dimensionsless  set: 


(2)  X  =  -g- 


Furthermore: 


Also  set: 


where 


Z  =  -T“ 


£ 


Ur  ttdn 

K  s  i  +  1  cos(e-cn) 


e>*f  =  zx 


U,=o 


Under  isothermal  conditions  \i  is  constant  and  the  dimensionless  Reynolds 
equation  becomes: 

W  £  [  ^  ix  ]  +  (' l)  di  [  $  $1  ]  =  12  n [&  +  (#  -  f )  sWe-«)] 

The  corresponding  finite  difference  equation  is: 


<5)f,r 


•V*  +  kVt  +(l)>(kJ[^/  kV  J 


where  the  finite  difference  mesh  is  as  follows: 


Figure  2. 
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The  pressure  et  the  bearing  end  is  ambient  (i.e.O).  At  the  leading  and 
the  trailing  edge  it  is  given  by  the  input  (partial  arc  only).  At  the 
bearing  centerline  there  is  symmetry.  Thus  the  boundary  conditions  to 
Eq.  (5)  become: 


(6) 


14  ^UelsM^0  >  tHth 


The  pressure  is  initially  set  equal  to  aero  and  Eq.  (S)  is  then  solved 


by  iterationa. 

The  load  carrying  capacity  is  calculated  from: 

ri  r** 


:r  *  -  2^  j  Pco$(e-d)  dxdz  *  -  ?  AX-AZ  ^  Pg  (oify-rf) 

2  M  Psin (e-g)  dxdz  -  2 ax •  ai £  Pg  sinfe^ 

JSu 


The  indicated  summations  are  computed  by  Simpson  rule  of  integration,  l.a. 

-  -  3  and  if  either  m  or  n  are  odd  the  flrat 

Interval  is  Integrated  by  T^A^P, +8lg~fg]  In  addition  end  corr¬ 
ections  are  made  at  the  boundary  to  the  ruptured  film. 


S-f  H  H  I  liW  V  lil-M  vTtTTw K 


The  convergence  of  the  pressure  iteration  la  tested  by  two  methods,  which 
shall  be  denoted  relative  convergence  and  absolute  convergence.  The  re¬ 
lative  convergence  is  tested  as  follows:  after  the  k'th  pressure  iteration 
compute: 


-?y- 


where  ds,  is  given  by  the  input.  When  Ag  becomes  zero  or  negative  then 
relative  convergence  has  been  achieved. 


To  compute  the  absolute  convergence  the  following 
After  each  iteration  the  sum  of  the  pressures  are 
suit  for  the  k'th  iteration  be  denoted: 

<“> 

Assume  that  after  infinitely  many  iterations  this 
and  set: 

%  -  -  A  e"  (a  w  B 

from  which: 

(11)  u  —  i ,  -4-  - 

Ljk  2 

is  calculated  after  each  Iteration  (exceptions,  see  later)  and  the 
absolute  convergence  is  computed  from: 

1  ~ 1  "  <S»  =  ^ 

where  6 ^  is  given  by  the  input.  When  becomes  zero  or  negative 

absolute  convergence  has  been  achieved. 


criteria  is  chosen, 
computed.  Let  the  re¬ 


sum  will  be  ijieo 
Constah(j) 


In  order  to  obtain  complete  convergence,  both  relative  and  absolute  con¬ 
vergence  must  be  satisfied.  Equations  (9),  (11)  and  (12)  are  printed  as 
output  after  each  iteration.  Eq.  (11)  serves  an  additional  purpose,  namely 
to  extrapolate  the  pressure  distribution.  When  becomes  a  smooth  curve 

and  starts  to  level  off  then  a  new  extrapolated  pressure  distribution  is 
calculated  from 


.  p“ 

*>J 


(13) 


and  the  pressure  iterations  proceed  from  this  new  distribution  either 
until  convergence  is  achieved  or  until  again  becomes  sufficiently 
smooth  that  Eq.(13)  may  be  used  once  more. 


It  remains  to  define  the  criteria  for  "smoothness".  As  applied  to  Eq.(ll) 
is  not  calculated  for  the  first.  2  iterations  or  for  the  2  itera¬ 
tions  following  a  pressure  extrapolation.  In  these  cases  the  computer 
output  shows: 


(14) 


=  !03fl 

Aa  -  1.0 


Furthermore,  4w~  is  not  calculated  if  the  -curve  is  not  monotoneiy 


increasing  with  decreasing  gradient.,  i.e.: 

V  4*-.  “4k-*  “4k)  >  0 

<15)  V  ( 1  *  '0'S.  J 

o.)  2  ,  (k- ^  2 

W  ^4k-i“4k-i'“4k)40 
(16)  c; 

<xy  (%-fabo  ] 

w  ( 2 - ** -4k)  >0 ,  lo\ j 


4koo  kij  e^(il) 

l/v38 

^ 

Aa  =  LO 
4k~»  “  4k 

&A  ”  o.o 


Similarly,  Eq.  (13)  is  not  applied  until  is  "smooth".  Basically 

this  is  determined  by  the  input  item  6&<  which  specifies  how  small  the 
relative  difference  between  two  consecutive  Iji^o  must  be  for  extrapola¬ 
tion  to  be  performed.  In  addition  the  -curve  must  be  increasing. 
Hence,  the  criteria  for  extrapolation  becomes: 
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(18) 


a)  0  *  *•  <*■« 

t;  0  *  ( )  4  ^ 

cj  A**  1.0 

k  *  -4  ,  ( k  “  ^«A^f*i4i*"^  ^ 


"\ 

> 


When  theee  equation* 
are  eetiefied  Eq»(13) 
is  performed,  otherwise 
not. 


These  Equations  (15)  to  (18)  ere  .11  used  In  the  704-ver.lon  of  the 
program.  In  the  1620-ver.ion  Eq.  (18,.)  1.  not  Included  end 
Eq.  (18, b)  1.  built  into  the  program,  given  the  v.lue  «»*-<»»  . 
more,  In  Eq.  (16)  Is  replaced  by  K> 


in 

Further 


Input 

The  input  form  is  shown  on  page  17  ■  Below  follow  the  instructions  for 
preparing  the  input. 

Card  1  and  2 

These  are  heading  cards  to  be  used  for  Identification  and  must  always  be 
given.  Punch  text  in  column  2  to  72  (2  to  52  in  1620  version). 

Surd  ,3 

L/P  Ratio  This  is  the  ratio  between  bearing  length  and  journal  diameter. 
Inlet  Angle  .  This  is  the  angle  in  degrees  measured  from  the  vertical 

reference  line  to  the  leading  edge  of  the  bearing  pad  in  the  direction  of 
Journal  rotation.  For  a  full  360°  bearing  ©  may  be  any  value  as  long 
as  “  360°. 

Outlet  Angle.  0nut  This  is  the  angle  in  degrees  measured  from  the 
vertical  reference  line  to  the  trailing  edge  of  the  bearing  pad  in  the  dir¬ 
ection  of  journal  rotation.  In  a  full  360°  bearing  @out  may  be  any 
value  as  long  as  Gout”®!*  *360°. 

Inlet  Pressure  Pin  This  is  the  dimensionless  pressure  along  the  lead¬ 
ing  edge  of  the  bearing  pad.  If  the  actual  pressure  1b  P  psig,  then 

Thus,  for  ambient  pressure  fjhs0  •  *n  t*ie  finite  difference  represen¬ 
tation  the  pressure  along  the  leading  edge  is  set  equal  to  f*,,  except 
at  the  end  of  the  bearing  (i»l)  where  the  pressure  is  ambient. 

Si 

fW.  This  is  the  dimensionless  pressure  along  the  trail¬ 
ing  edge  of  the  bearing  pad.  For  consents,  see  above. 

RglgUvi  CPtivergance  Limit  6a  This  limit  determines  when  relative  con¬ 
vergence  has  been  achieved  i.e.  when  the  relative  difference  between  two  con¬ 
secutive  pressure 'iterations  is  less  than  Jg  .  For  details,  see  Eq.  (9) 
page  6  ,  A. typical  value  is  ($^*,00 1  or  .002 


Absolute  Convergence  Limit  oa  This  limit  determines  when  absolute 
convergence  has  been  achieved,  i.e.  when  the  difference  between  a  given 
pressure  iteration  and  the  corresponding  extrapolated  final  value  be¬ 
comes  less  than  6 \  .  For  details,  see  Eq.  (12),  page  6*  A  typical 

value  may  be  (fAs.0O2  to  .005.  It  should  be  noted  that  absolute  con¬ 
vergence  is  much  more  strict  than  relative  convergence  i.e.  6^  should 
be  2  or  3  times  greater  than  for  consistency.  Hence,  6^=  ,002 

Is  rather  small,  especially  if  the  mesh  has  many  points  (150  to  200  and 
more).  Considering  that  <f^-,002  means  an  absolute  accuracy  of  .2  per¬ 
cent  which  in  most  cases  is  much  finer  than  needed,  and  also  that  the 
method  of  solution  is  approximate  anyway,  6 should  not  be  made  any 
smaller  than  actually  needed.  The  number  of  pressure  iterations  increases 
rapidly  with  decreasing  6 \  and  for  most  practical  cases  <^  =  .00?  should 
give  sufficient  accuracy. 

Relaxation  Factor  for  Pressure  Calculation  In  order  to  increase 

(  or  decrease)  the  rate  of  convergence  of  the  pressure  iterations  a  re¬ 
laxation  factor  may  be  used.  To  illustrate  assume  that  in  the  k’th  itera- 
D(k) 

tion  the  pressure  rj ,  has  been  computed  from  Eq.  (5)  for  any  meshpoint. 

/;  0W 

Then  the  actually  stored  value  is  not  r^j  but 


Hence,  will  accelerate  the  convergence  and  will  decelerate 

the  convergence.  Setting  fp® I  results  in  no  relaxation.  It  should  be 
noted  that  the  program  includes  an  additional  feature  for  speeding  up  the 
calculation,  namely  the  pressure  extrapolation  given  by  Eq.  (13).  If  ex¬ 
trapolation  is  desired  it  is  questionable  if  relaxation  should  also  be  used. 


Card  5 


1.  Number  of  Circumferential  Subdivisions ,m  As  shown  in  Figure  2,  the 
circumferential  length  of  the  bearing  is  subdivided  into  m  increments 

aw  —  -ZL-  lSemjz.Sin.1  , 

AX  “  180  *rt  •  Thus  the  number  of  meshpoints  in  the  circumferential 

direction  is  (nrt-1) . 

2.  Number  of  Axial  Subdivisions,  n  As  shown  on  Figure  2,  the  bearing  half- 

length  is  subdivided  into  n  increments  .  Thus  the  number  of  mesh¬ 

points  in, the  axial  direction  is  (n+l). 

3.  Full  or  Partial  Bearing  If  this  item  is  0  the  program  assumes  that  the 
bearing  is  full  360°  such  that  the  pressure  along  the  leading  edge  is  equal 


to  the  pressure  along  the  trailing  edge  with  a  continuous  gradient.  Hence, 
item  4,  Card  3,  and  item  l,Card  4,  are  ignored. 

If  this  item  is  1  the  program  assumes  that  the  bearing  is  partial  arc. 

4.  Humber  of  Pressure  Iterations  In  order  to  ensure  that  the  computer 
does  not  use  an  excessive  amount  of  time  to  meet  the  convergence  limits 
this  item  sets  a  limit  for  the  maximum  number  of  pressure  Iterations. 

If  this  number  is  exceeded  the  computer  gives  the  results  obtained  up 
to  this  point. 

5.  Humber  of  £-tf- Input  Frequently  the  load  and  other  quantities  are 
wanted  as  a  function  of  eccentricity  ratio  E,  keeping  other  parameters  such 
as  L/D  ratio  and  arc  length  the  same.  The  program  provides  for  giving  as 
many  values  of  t  as  desire!  and  the  present  item  specifies  how  many  cases 
there  are,  i.e.  how  many  of  Card  6  (or  of  pairs  of  Card  6-Card  7). 

6.  Input  Control  If  this  item  is  zero  the  program  will  return  after  com¬ 

pleted  calculation  to  read  in  an  additional  input  set.  If  this  item  is  1 
the  program  assumes  this  to  be  the  last  input  set  and  will  go  to  normal 
stop  after  completed  calculation.  Hote,  that  this  of  course  does  not  refer 
to  the  above  mentioned  E-A- input  which  is  regarded  as  a  part  of  an  input 
set.  An  input  set  starts  with  Card  1. 

7.  Extrapolation  Limit.  <fsx  (704-version  only).  After  each  pressure 

iteration,  say  the  k'th,  all  pressures  are  summed  up.  Denote  the  sum  . 

Based  on  3  consecutive  values  of  an  exponential  extrapolation  is 
used  to  calculate  the  value  of  the  sum  after  infinitely  many  Iterations, 
denoted  Ijk*,  .  When  the  relative  difference  between  three  consecutive 
values  of  is  less  than  an  extrapolation  of  the  pressure  dis¬ 

tribution  is  performed  by  multiplying  each  pressure  by 

details,  see  Eq.  (13)  and  (18),  page  6«4J.  A  few  remarks  are  needed.  It 
is  obvious  that  the  intention  of  the  extrapolation  is  to  speed  up  the  con¬ 
vergence  of  the  pressure  iteration  but  it  is  difficult  to  ascertain  that 
this  is  always  the  case.  When  the  absolute  convergence  limit  6^  is  not 
too  tight  the  extrapolation  may  reduce  the  number  of  iterations  by  more 
then  30  percent  but  when  <fA  is  small  it  is  more  difficult  to  evaluate 
the  effect.  After  an  extrapolation  the  pressure  distribution  is  no  longer 
in  "equilibrium"  and  this  reflects  clearly  in  the  pressure  sums  and 


very  strongly  in  the  extrapolated  sums  In  general  \jk  will  be 

decreasing  right  after  an  extrapolation  such  that  cannot  be  computed 

and  therefore  is  given  the  arbitrary  value  10  (10  in  1620-verson)  in 

the  output  .  After  a  certain  number  of  iterations  tj*  will  slowly  start  to 
increase  again  and,  when  smooth  enough,  a  new  extrapolation  will  be  perform¬ 
ed,  repeating  the  cycle.  When  6 \  is  small  and  the  number  of  meshpoints  is 
large  these  cycles  may  be  rather  long  which  tends  to  reduce  the  benefit  of 
the  extrapolation.  If  it  is  not  desired  to  extrapolate,  set  If 

extrapolation  is  desired,  an  example  of  its  value  is  <Se x  =  .005. 


Card  6 


There  should  be  as  many  Card  6  as  given  by  item  5,  Card  5.  Each  Card  6 
with  MATT=+1  (item  5  below)  should  be  followed  by  a  Card  7. 

Eccentricity  Ratio.  £  This  is  the  ratio  between  the  journal  eccentri¬ 
city  and  the  radial  clearance,  see  Figure  1. 

Attitude  Angle,  Ot  This  is  the  angle  in  degrees  measured  from  the 
vertical  reference  line  to  the  line  of  centers  (i.e.  the  line  connecting 
bearing  and  journal  centers)  in  the  direction  of  rotation,  see  Figure  1. 

Radial  Squeeze  Film  Velocity.  ^CO _ This  is  the  dimensionless  velocity 

of  the  journal  center  along  the  line  of  centers,  see  Figure  1.  However,  if 
item  5,  Card  6,  is  MATI--1  then  this  item  means  for  use  in  calcula¬ 

ting  force  derivatives  with  respect  to  ^ 

Tangential  Squeeze  Film  Velocity,  £<Va>  This  is  the  dimensionless  velocity 
of  the  journal  center  perpendicular  to  the  line  of  centers,  see  Figure  1. 

MATT  If  MATT=-1  a  Card  7  must  be  given.  The  program  will  perform  a  number 
of  complete  pressure  calculations,  changing  the  attitude  angle  in  order  to 
make  the  horizontal  force  component  zero  .  Thereafter  an  additional  4  complete 
calculations  will  be  performed  to  calculate  and 

If  Matt-!0  no  Card  7  is  required.  Only  one  complete  calculation  will  be  per¬ 
formed  based  on  the  input  value  for  the  attitude  angle. 

If  MATT=+1  a  Card  7  is  required  after  Card  6.  The  program  will  iterate  on  the 
attitude  angle  in  order  to  eliminate  the  horizontal  force  component.  However, 
force  derivatives  will  not  be  computed. 
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Card  7 

Card  7  can  only  be  given  when  MATT*+1  (see  item  5  above) 

Eccentricity  Ratio  Increment  AS  When  MATO+l  this  item  is  not  used. 

When  MATT®-1  A£  is  used  in  computing  the  derivatives  with  respect  to  t. 

Upon  completion  of  the  attitude  angle  iterations  the  program  proceeds  to 
calculate  the  change  in  bearing  load  as  a  function  of  £  ,  keeping  all  other 
parameters  constant.  Denoting  the  original  eccentricity  ratio  £0  (item  1, 

Card  6)  the  program  performs  calculations  for  £S£0+AE  and  for  £“£,**A£. 

Hence 3  the  results  may  be  used  to  calculate  64/dt  = 

Convergence  Limit  for  Attitude  Angle  Iteration,  .  <5^  is  in  degrees  and 
determines  when  the  attitude  angle  iteration  has  converged.  After  the  k'th 
iteration  the  program  computes  the  value  of  the  attitude  angle  as 

t  t  *  W(VO 

If  | ci<k+,)  “  otlh'  |  —  6+  the  calculation  has  converged  (or  if  item  4  below  is 
exceeded)  and  the  program  proceeds  to  calculate  derivatives.  Otherwise  a 
new  trial  value  of  tk  is  determined  as  described  below  and  a  new  calculation 
is  performed.  A  typical  value  is 

Relaxation  Factor  for  Attitude  Angle  Iteration  As  described  above/ a 

new  pressure  calculation  is  performed  when  k'-'W’K .  This  new  calcula¬ 
tion  is  based  on  an  <A-value  determined  as  follows:  the  first  calculation 
uses  the  original  value  for  ol  (item  2,  Card  6).  The  two  next  calculations  use 

Any  additional  calculations  uses  an  ci  -value  determined  by  parabolic  inter¬ 
polation.  ■£  is  actually  a  function  of  eccentricity  ratio,  L/D-ratio,  bear¬ 
ing  geometry  etc.  It  usually  has  a  maximum  around  £  =.5  where  £=2.5  to  3. 

In  general,  use  ■£  =U),  to  1.5. 

Maximum  Number  of  Attitude  Angle  Iterations  This  item  limits  the  number  of 
attitude  angle  iterations  so  that  an  excessive  number  of  iterations  is  not  per¬ 
formed  even  if  convergence  has  not  beei  obtained.  Normally  5  to  6  iterations 
are  needed  if  <£  *.01°. 
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Outcut 


The  output  is  in  general  self-explanatory.  First  come  the  two  heading  cards 
followed  by  3  or  4  lines  listing  the  input.  Thereafter  the  convergence  of  the 
pressure  iteration  is  printed  as  a  list  with  6  columns.  The  first  column 
gives  the  number  of  the  iteration,  the  second  column  is  Ag  from  Eq.  (9), 

the  third  column  is  Aa  from  Eq.  (12),  the  fourth  column  is  ^  from  Eq.(10), 

the  fifth  column  is  the  extrapolated  sum  y* from  Eq.  (11)  or  from  Eq.  (15) 

to  (17)  and  the  sixth  column  gives  a  count  of  the  number  of  iterations  after 

last  pressure  extrapolation.  In  the  1620-version  column  5  and  6  are  inter¬ 
changed.  Note,  that  for  convergence  Ag  and  A^  must  be  zero  or  negative. 

Next  follows  the  pressure  distribution  listed  as: 


rota-tion 


IU 


— ‘mim _ up _ 

- 1 

1 

1 

_ trajlm _ edit _ 

n 

i  .  , 

_ edit _ 1 

1  '  1 

1 

J 

1 

k 

1 

1 

1 _ ’traa«»|  tdje 

r 

• 

The  output  format  provides  for  a  maximum  of  7  columns.  Hence,  if  (n+1)  ex¬ 
ceeds  7  the  pressure  field  is  broken  up  as  shown  above.  The  last  line  gives 
the  dimensionless  radial  force  component  fr  (Eq.(7)),  the  tangential  force 
component  (Eq.(8)),  the  vertical  force  component  fv  ,  the  horizontal  force 


the  total  force 

where: 


j  the  Sommerfeld  NumberS  and  the  cal- 


component  ^ 

culated  attitude  angle^ot^i, 

fv  *  J^Pcose  dx  dz  -  fr&s«i  + 

^  f^smedxdz  =•  — fhSm<A 

*  ft?  +  C 

c  '  :  jdL  (JLf 

^  f/dl  V  c  / 

<*c.k  =  W  (Vfr) 


Hence,  to  convert  the  dlnensionleas  force  f  to  ectuel  force  7; 

F»^lDLN(?)f  Iks 

Similarly,  the  actual  preasure  F  in  psig  is  obtained  from  the  dimensionless 
pressure  7  by 

P  ”  /a  N  (£)’  P  pstg 

Computer  Operation 

The  input  is  on  cards.  In  the  704-version  the  output  is  given  on  Tape  3, 
vhich  is  the  only  tape  used  by  the  program.  Upon  completion  of  calculation 
Tape  3  is  given  an  "End  of  File",  but  not  rewound,  and  the  computer  stops 
with  77777g  in  the  address  field.  The  stop  is  at  3361g.  The  program  is 
compiled  with  part  of  LSTG's  TOF-system,  i.e.  senseswitch  1  and  2  down. 
The  1620-version  is  written  in  FORTRAN  I,  both  input  and  output  is  on 
cards. 


Computer  Time 


The  two  most  important  factors  in  determining  the  computer  time  is  the 

number  of  meshpoints  m  x  n  and  the  absolute  convergence  limit  6: \  . 

Whereas-,  it  is  difficult  to  evaluate  the  effect  of  the  latter,  it  is 

approximately  correct  to  state  that  the  computer  time  is  proportional  to 
2 

(m  x  n)  .  This  is  due  to  two  contributions:  a)  for  each  iteration  the 
flme  is  of  course  directly  proportional  to  the  number  of  meshpoints.  This 
time  may  be  estimated  from: 

1000  mesh  point  calculations  (i.e. of  )  per  5  seconds, 
b)  the  number  of  Iterations  needed  for  convergence  increase*,  approximately 
proportional  to  the  number  of  meshpoints.  Example:  for  6 a  *  .002  and 
(SgX  -  .005  the  average  number  of  iterations  is: 


No.  iterations 
m  x  n  ■  12  x  6  46.6 

m  x  n  ■  16  x  6  65 

m  x  n  ■  36  x  10  Approx. 200 


time _ 

17  seconds 
31  seconds 
6  minutes 
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The  first  two  numbers  sre  averaged  over  several  values  of  e  .  It  should 
be  noted  that  €  also  affects  the  number  of  Iterations.  As  an  example, 
take  the  same  conditions  as  above.  Then: 

mxn«12x6  n  x  n  <  16  x  6 

No.  iterations  No  iterations 


e 

a 

.1 

44 

70 

€ 

m 

.2 

48 

75 

e 

m 

.4 

54 

74 

E 

m 

.6 

47 

55 

E 

m 

.8 

31 

55 

From  the  above  both  computer  time  and  number  of  iterations  may  be 
estimated. 


INPUT  FOR  PN0091  (IBM  704) 
HYDRODYNAMIC  INCOMPRESSIBLE  JOURNAL  BEARING 


Card  1.  Text,  Col.  2-72: 


Card  2.  Text.  Col.  2-72: 


Card  3 
(1P4E15.7) 

E _ _  L/D  -Ratio 

E  _  0in>  Inlet  Angle, degrees 

E  _  ®out*  ®ut*et  Angle , degrees 

E  _  ^in*  *n^et  Pressure,  dimensionless 


Card  4 
(1P4E15.7) 

E  P  ,  Outlet  Pressure,  dimensionless 

E _ R,  Relative  Convergence  Limit 

E _ <$A,  Absolute  Convergence  Limit 

E  _  fp,  Relaxation  Factor  for  Pressure 

Calculation. 


Card  5 

(615 , 1PE15. 7) 

1.  m.  Number  of  Circumferential  Subdivision  (m  4  200) 

2.  n,  Number  of  Axial  Subdivision  on  Half  length  (n 100) 

3.  0:  Full  360°  Bearing  1:  Partial  Arc  Bearing 

4.  Maximum  Number  of  Pressure  Iterations 

5.  Number  of  £-ct-  Input 

6.  0:  More  Input  sets  follow  1:  Last  Input  set 

_ _ E _ 7.  /  Extrapolation  Limit 
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Card  6 

(1P4E15.7,  15) 

_ E _ £  ,  Eccentricity  Ratio 

_ E _ oC  ,  Attitude  Angle,  degrees 

_  E  _  Radial  Squeeze  Film  Velocity , dimensionless 

_  E  _  Tangential  Squeeze  Film  Velocity , dimensionless 

MATT,  *1:  Attitude  Angle  Iteration  4-  Force  Derivatives 
0:  Single  Calculation 

+1:  Attitude  Angle  Iteration,  no  Derivatives 


Card  7 

(1P3E15.7,  15) 

Only  if  Card  6,  item  5,  MATT-+1 


E  _  >  Eccentricity  Ratio  Increment  for 


t 


E  _  <jj  ,  Convergence  Limit  for  Attitude  Angle  Itera- 

tion,  degrees. 

E  ■£  ,  Relaxation  Factor  for  Attitude  Angle  Itera- 

*  tion 

_  Maximum  Number  of  Attitude  Angle  Iterations. 


Note:  An  Input  Set  consists  of  Card  1  to  Card  5  plus  as  many  pairs  of  Card  6- 
Card  7  as  given  by  item  5,  Card  5. 


INPUT  FOR  PN0091  (IBM  1620) 
HYDRODYNAMIC  INCOMPRESSIBLE  JOURNAL  BEARING 
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Car<jt  1.  Text,  Col.  2-52: 


Card  2.  Text,  Col.  2-52: 


Card  3 

(1XE14.7,  1XE14.7,  1XE14.7,  1XE14.7) 

_ _ _ E _ L/D  -Ratio 

_ _ E _ ®in*  *n*,et  Angle,  degrees 

_ _ E _ ®out*  Angle,  degrees 

_ _  E  Pjn>  Inlet  Pressure,  dimensionless 

Card  4 

(1XE14.7,  1XE14.7,  1XE14.7,  1XE14.7) 

F.  P  .  ,  Outlet  Pressure,  dimensionless 

■ - -  out 

_ _ E _ Relative  Convergence  Limit 

_ _ E _ Absolute  Convergence  Limit 

_ _ E _ Relaxation  Factor  for  Pressure 

Calculation. 

Card  5 

(1X14,  1X14,  1X14,  1X14,  1X14,  1X14) 

1.  m.  Number  of  Circumferential  Subdivision  (m  4-  26) 

2.  n.  Number  of  Axial  Subdivision  on  Half  length  (n  6  12) 

3.  0:  Full  360°  Bearing  1:  Partial  Arc  Bearing 

4.  Maximum  Number  of  Pressure  Iterations 

5.  Number  of  £-*- Input 

6.  0:  More  Input  Sets  follow  1:  Last  Input  Set 
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Sft&iL 

(1XE14.7,  1XE14.7,  1XE14.7,  1XE14.7,  1X14) 

_ _ E _  £  ,  Eccentricity  Ratio 

_ _ —  *  si,' Attitude  Angle,  degrees 

_  _ _ _  e _ „  _  Radial  Squaese  Film  Velocity,  dimensionless 

_ _ _ _ _  _  e  -  —  -  ^5  »  Tangential  Squaese  Film  Velocity,  " 

_  MATT,  *1:  Attitude  Angel  Iteration  +  Force  Deri- 
0:  Single  Calculation.  vatives. 

'•+1:  Attitude  Angle  Iteration,  no  Derive  - 

tlves. 


for.4.7 

(1XB14.7,  1XE14.7,  1XE14.7,  1X14) 

Only  if  Card  6,  Item  5,  Matt-+1 

M 

_ g _ A£»  Eccentricity  Ratio  Increment  for 

E  <L  ,  Convergence- Limit  for  Attitude  Angle  Itera- 

”  ~  ,  tlon.  degrees 

_  E  %  •  Relaxation  Factor  for  Attitude  Angel  Itera- 


_  Maximum  Humber  of  Attitude  Angle  Iterations 


Rote:  An  Input  Set  consists  of  Card  1  to  Card  5  plus  as  many  pairs  of 
<Jard  6*  Card  7  as  given  by  item  5,  Card  5. 
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TABU  I 


D  B  ®  *  a  fr 

£ 

5(A) 

'01 

& 

a  (A) 

'Oi 

25°ABC  1640.0  .01  84.79  5.55xl0'3 

.00558  .00161 

6.09x10 

'4  .0610 

.0111 

-.0601  .00937 

159.0  .1  65.00  .00266 

.0302 

.0257 

.00571 

.0645 

.0533 

-.0514  .0370 

70.7  .2  50.81  .00894 

.0641 

.0744 

.011 

.0781 

.0893 

-.0585  .0414 

.25  1.15  29.5  .35  37.04  .027 

.167 

.208 

.0204 

.125 

.154 

-.0198  .103 

12.5  .5  28.17  .0706 

.470 

.542 

.0378 

.249 

.282 

.0226  .183 

4.53  .65  21.22  .206 

1.75 

1.70 

.0799 

.667 

.633 

.187  .401 

1.06  .8  15.03  .906 

12.3 

9.22 

.243 

3.22 

2.28 

1.32  1.36 

.0459.95  8.34  21.5 

1000.0 

428.0 

3.16 

134.0 

44.9 

64.7  19.9 

50°  ABC  476.0  .01  87.38  9.55xl0*S 

.00979 

.0131 

.0021 

.21 

.0177 

-.203  .013 

45.8  .1  71.23  .00703 

.0799 

.06  33 

.0207 

.228 

.14 

-.194  .0813 

20.4  .2  56.84  .0268 

.186 

.2 

.0410 

.275 

.268 

-.171  .155 

.25  .573  8.63  .35  42.78  .0851 

.502 

.593 

.0787 

.445 

.487 

-.146  .272 

3.76  .5  33.29  .222 

1.37 

1.56 

.146 

.85 

.886 

-.0874  .472 

1.44  .65  25.88  .625 

4.87 

4.92 

.303 

2.2 

1.94 

.19  .934 

.383  .8  19.68  2.46 

29.2 

24.8 

.879 

9.46 

6.0 

2.31  2.37 

.0262.95  12.15  37.3 

1410.0 

611.0 

8.04 

247.0 

76.0 

27.3  12.3 

209.0  .01  86.76  2.69xl0'4 

.0274 

.0219 

.00478 

.479 

.0532 

-.467  .0402 

20.3  .1  67.19  .0191 

.217 

.185 

.0454 

.505 

.381 

-.413  .256 

9.16  .2  52.14  .067 

.476 

.565 

.0862 

.597 

.669 

-.317  .443 

3.91  .35  38.44  .2 

1.2 

1.57 

.159 

.919 

1.14 

-.139  .741 

3.0  .4  35.18  .272 

1.61 

2.12 

.192 

1.09 

1.37 

-.0448  .89 

.5  1.15  1.73  .5  29.68  .501 

3.17 

3.93 

.286 

1.73 

2.0 

.165  1.22 

.951  .6  25.13  .952 

6.76 

7.81 

.447 

3.0 

3.17 

.707  1.88 

.678  .65  23.11  1.36 

10.6 

11.6 

.579 

4.26 

4.18 

1.26  2.38 

.188  .8  17.63  5.06 

59.1 

53.9 

1.61 

17.2 

12.5 

7.48  6.02 

.0158  .95  11.79  62.0 

2111.0 

1040.0 

12.9 

365.0 

129.0 

68.8  27.7 

157.0  .01  86.16  4.25x10' 

4  .0431 

.0271 

.00636 

.637 

.0845 

-.6  23  .  0667 

15.3  .1  65.0  .0276 

.315 

.287 

.0591 

.663 

.551 

-.518  .388 

6.97  .2  49.92  .0924 

.662 

.831 

.11 

.77 

.923 

-.359  .643 

75  1.72  3.0  .35  36.65  .267 

1.61 

2.21 

.199 

1.16 

1.53 

-.079  1.03 

1.34  .5  28.25  .656 

4.17 

5.4 

.353 

2.15 

2.62 

.459  1.7 

.53  .65  22.03  1.75 

13.6 

15.6 

.708 

5.2 

5.38 

2.17  3.24 

.15  .8  16.90  6.36 

72.5 

70.1 

1.93 

20.3  15.9  11.0  8.08 

.0136  .95  11.7  71.9 

2360.0 

L  230.0 

14.9 

407.0  150.0  91.9  34.6 

60°  ARC  371.0  .01  87.85  1.01x10 

**  .0105 

.0248 

.00269 

.269 

.0167 

-.2*  .0124 

35.8  .1  72.11  .00859  .0981 

.0876 

.0266 

.289 

.173 

-.246  .0923 

16.0  .2  58.45  .0327 

.228 

.25 

.0533 

.352 

.324 

-.224  .178 

.25  .477  6.79  .35  44.54  .105 

.612 

.740 

.103 

.564 

.603 

-.191  .313 

3.01  .5  35.0  .272 

1.66 

1.94 

.19 

1.07 

1.08 

-.11  .533 

1.19  .65  27.76  .746 

5.71 

5.98 

.393 

2.7 

2.28 

.224  .96 

.333  .8  21.67  2.79 

31.6 

28.8 

1.11 

11.5 

6.94 

2.18  2.17 

.0249  .95  13. 11  39.1 

1440.0 

620.0 

9.11 

266.0 

81.7 

20.0  9.39 

TABLE  I  (Cont.) 


1 

D 

1 

B 

S  c 

a 

fr 

£ 

&<*) 

0) 

£ 

s<*> 

CD 

£ 

60° 

ARC 

147.0  .01 

87.33 

3.18xl0*4 

.0327 

.0452 

.00682 

.683  .0597 

-.659 

.044 

(Cone'd) 

14.2  .1 

68.43 

.0259 

.295 

.261 

.0655 

.73 

.52 

-.591 

.336 

6.42  .2 

53.9 

.0918 

.647 

.764 

.126 

.854  .917 

-.476 

.584 

.5 

.955 

2.77  .35 

40,26 

.276 

1.6 

2.12 

.233 

1.29 

1.57 

-.266 

.963 

1.25  .5 

31.38 

.684 

4.21 

5.3 

.417 

2.41 

2.73 

.144 

1.58 

.502  .65 

24.84 

1.81 

13.4 

15.4 

.836 

5.7 

5.54 

1.43 

2.85 

.149  .8 

19.54 

6.34 

68.8 

68.5 

2.25 

21.8 

15.9 

8.6 

6.32 

.0143  .95 

12.78 

68.0  2150.0  1110.0 

15.4 

394.0 

141.0 

47.0 

23.1 

103.0  .01 

87.17 

4.78x10 

*4  .049 

.0543 

.00968 

.969 

.0935 

-.941 

.0686 

10.1  .1 

66.49 

.0396 

.452 

.411 

.091 

1.01 

.792 

-.804 

.545 

4.6  .2 

51.32 

.136 

.967 

1.05 

.17 

1.17 

1.36 

-.573 

.922 

.75 

1.43  2.0  .35 

38.09 

.393 

2.32 

3.2 

.308 

1.74 

2.24 

-.182 

1.46 

.913  .5 

29.69 

.951 

5.85 

7.75 

.542 

3.15 

3.8 

.547 

2.34 

.373  .65 

23.55 

2.45 

18.2 

21.9 

1.07 

7.33 

7.6 

2.72 

4.18 

.115  .8 

18.34 

8.25 

87.3 

93.2 

2.77 

26.6 

20.4 

13.7 

9.1 

.0121  .95 

12.68 

80.3 

2420.0  : 

1320. 0 

18.1  440.0 

166.0 

63.5  30.1 

Ml 

ARC  90.7  .01 

87.91 

4.02x10* 

4  .0429 

.137 

.011 

1.1 

.0621 

-1.03 

.05 

8.79  .1 

71.31 

.0365 

.424 

.431 

.108 

1.17 

.729 

-.951 

.435 

3.99  .2 

57.35 

.135 

.95 

1.16 

.211 

1.37 

1.35 

-.796 

.788 

.5 

.716  1.76  .35 

43.94 

.41 

2.34 

3.17 

.395 

2.05 

2.34 

-.512 

1.28 

.821  .5 

35.14 

.995 

5.76 

7.85 

.701 

3.61 

3.95 

.0624 

1.94 

.351  .65 

28.81 

2.49 

17.4 

22.0 

1.37 

8.31 

7.52 

1.82 

2.99 

.116  .8 

23.33 

7.94 

78.8 

80.2 

3.42 

28.4 

20.6 

5.86 

5.47 

.013  .95 

14.27 

74.5 

2220.0 

1140.0 

19.0 

433.0  154.0 

24.4 

15.0 

57.4  .01 

87.73 

6.92xl0*4 

.073 

.179 

.0174  1.74 

.114 

-1.65 

.0931 

5.61  .1 

68.76 

.0645 

.745 

.738 

.166 

1.81 

1.29 

-1.45 

.849 

2.57  .2 

54.47 

.226 

1.59 

2.01 

.316 

2.08 

2.26 

-1.1 

1.44 

.75 

1.07  1.15  .35 

41.37 

.654 

3.71 

5.28 

.576 

3.01 

3.73 

•  .488 

2.24 

".  901  .4 

38.17 

.873 

4.85 

7.05 

.686 

3.49 

4.39 

-  .211 

2.56 

.546  .5 

33.0 

1.54 

8.87 

12.5 

.998 

5.21 

6.13 

.622 

3.34 

.321  .6 

28.89 

2.73 

17.2 

23.4 

1.5 

8.4 

8.96 

2.2 

4.46 

.24  .65 

27.27 

3.71 

24.9 

33.4 

1.92 

11.1 

11.3 

3.5 

5.19 

.0832  .8 

22.50 

11.1 

103.0 

116.0 

4.6 

36.4 

27.7 

11.2 

8.87 

.0108  .95 

14.29 

89.8 

2480.0 

1370.0 

22.9 

492.0  186.0 

31.6 

21.2 

45.5  .01 

87.46 

9.73x10* 

4  .102 

.203 

.022 

2.2 

.171 

-2.09 

.155 

4.46  .1 

67.49 

.0859 

.989 

.979 

.207 

2.27 

1.72 

-1.78 

1.18 

2.06  .2 

52.84 

.293 

2.06 

2.69 

.387 

2.56 

2.93 

-1.25 

1.96 

1.0 

1.43  .931  .35 

39.89 

.824 

4.67 

6.86 

.689 

3.61 

4.71 

-  .329 

2.96 

.488  .5 

31.90 

1.9 

10.8 

15.8 

1.18 

6.13 

7.58 

1.15 

4.33 

.20  .65 

26.33 

4.48 

29.5 

41.1 

2.22 

13.2 

13.7 

4.98 

6.66 

.0712  .8 

21.98 

13.0 

118.0 

139.0 

5.25 

41.3 

32.4 

15.6 

11.1 

.00992  .95  14.29 

97.7 

2640.0 

1490.0 

24.9 

514.0  202.0 

36.2 

24.7 

CRIT.  MASS 


Fig.  1  Six  50°  Tilting  Pads,  Centrally  pivoted,  L/D  -  25 
L/B  -  .573,  C'/C  -  1.  Load  between  pad*.  No  pad 
lntertia. 


Six  50°  Tilting  P«d»,  Centrally  pivoted,  L/D  -  , 
L/B  -  1.146,  c’/c  -  1.  Load  Mo  Pad 


Inertia . 


Fig.  3  Six  50°  Tilting  Pads, 
L/B  -  1.719,C'/C  -  1. 
Inertia. 


Centrally  pivoted,  L/D  -  .7! 
Load  between  pads.  No  pad 


CwC 


-79- 


Fig.  4  Five  60°  Tilting  Pads, 
L/B  =  .477,  C’/C  =  1. 

inertia. 


Centrally  pivoted,  L/D  =  .25, 
Load  between  pads.  No  pad 


CWMcrjt 
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Fia.  5  Five  60°  Tilting  Pads, 
L/B  -  .955,  C'/C  -  1. 
inertia. 


Centrally  pivoted,  L/D  ■  .5, 
Load  between  pads.  Ho  pad 


s=  w  ic) 

Four  80°  Tilting  Pads,  Centrally  pivoted,  L/D  -  .7 
L/B  -  1.074,  C'/C  -  1.  Load  between  pad*.  No  pad 
Inertia. 
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Fads,  Centrally  pivoted,  L/D  “  1.0, 
;  .  1.  Load  between  pads.  No  pad 


CWMcrit 
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Fig.  13  Four  80  Tilting  Pads, 
L/B  -  1.074,  C'/C  -  1. 
of  pad  inertia. 


Centrally  pivoted,  L/D  *  .75, 
Load  between  pads.  Effect 
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Fig.  14 


Four  80°  Tilting  Pads, 
L/B  -  1.074,  C'/C  -  1. 


Centrally  pivoted,  L/D  ■  .75, 
Load  betveen  pads.  Effect  of 
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Fig.  15  Four  80°  Tilting  Pads,  Centrally  pivoted,  L/D  -  .75, 
L/B  ■  1.074,  C'/C  ■  1.  Load  between  pads.  Effect 
of  pad  inertia. 
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Fig.  16  Four  80°  Tilting  P«ds,  Centrally  pivoted,  L/D  -  .75, 
L/B  -  1.074,  C'/C  -  1.  Load  between  pads.  Effect 
of  pad  Inertia. 


-92- 


C 

Fig.  17  Twelve  25°  Tilting  Pads,  Centrally  pivoted,  L/D  ■  .25, 

L/B  «  1.146.  Vertical  rotor.  Effect  of  preload  (l-C'/C). 


NMcrit 


Ftg.  21  Four  80°  Tilting  Pads,  Centrally  pivoted,  L/D  ■  .75, 

L/B  ■  1.074,  C'/C  «  1.  Load  between  pads.  Couparlson 
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NOMENCLATURE 


c 

c' 


CxxjCx^G^Cfl 

D 


e 

*• 


F 

f 

f?*S 

I 

KW/Kx^K,X/K„ 

L 

M 


Met 

N 


P4 

R 


Pad  clearance  (radius  of  curvature  of  pad  minus  journal  radius)inch. 

Pivot  circle  clearance  (radius  of  pivot  circle  minus  journal  radius) 

inch 

Bearing  damping  coefficients,  lbs. sec/in. 

Fixed  pad  damping  coefficients,  lbs. sec/in. 

Tilting  pad  damping  coefficients,  lbs. sec/in. 

Journal  diameter,  inch 

Journal  center  eccentricity  with  respect  to  pad  center,  inch 
Journal  center  eccentricity  with  respect  to  bearing  center,  inch 
Load  on  pad.  lbs. 

Radial  and  tangential  components  of  pad  load,  lbs. 

Components  in  £  and  /^-direction  of  pad  load, lbs (see  Fig.  22) 

=  F/Acj,  dimensionless  pad  force 

Radial  and  tangential  components  of  dimensionless  pad  force 
Components  in  f  and  /^-direction  of  dimensionless  pad  force 

2 

Transverse  mass  moment  of  inertia  of  shoe  around  pivot,  lbs. in. sec 

Bearing  spring  coefficients,  lbs/in. 

Fixed  pad  spring  coefficients,  lbs/in. 

Tilting  pad  spring  coefficients,  lbs/in. 

Bearing  length,  inch 
2 

=  l/R  p,  equivalent  pad  mass,  lbs. sec*/ in 

Value  of  equivalent  pad  mass  to  cause  pad  motion  resonance,  lbs.secVin 
Rotational  speed  of  journal,  RPS 
Coefficients  defined  by  Eq.  (31)  and  (32) 

Journal  radius ,  inch 
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S 

w 

£ 

£o 

> 

X 

/* 

<p 

V 

(Jj 


Radius  from  pad  center  to  actual  pivot  point  of  pad,  inch 
2 

*  (uNDL/W)  •  (R/C)  j  bearing  Sommerfeld  number 
2 

=  (jiNDL/F)  •  (R/C)  pad  Sommerfeld  number 
Bearing  load,  lbs . 

Coordinates  of  journal  center  with  respect  to  the  bearing. 

See  Fig.  22,  inch 

=  e/C,  eccentricity  ratio  with  respect  to  the  pad  center 

=  eQ/C,  eccentricity  ratio  with  respect  to  the  bearing  center 

Amplitude  for  pad  center  motion,  see  Fig.  22,  inch 

Amplitude  of  the  center  of  a  massless  pad,  inch. 

=  (nRL/ft )■ (R/C)^  bearing  coefficient 

2 

Lubricant  viscosity,  lbs. sec. /in 

Coordinates  of  journal  center  with  respect  to  the  pad,  see  Fig. 22 

inch 

Attitude  angle  with  respect  to  the  pad  load  line,  radians 

Attitude  angle  with  respect  to  the  bearing  load  line,  radians 

Angle  from  vertical  (negative  x-axis)  to  pad  pivot  point,  see  Fig.  22 

inch 

Angular  speed  of  shaft,  radian/sec. 


